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Abstract 



We investigate the asymptotic behavior of a generahzed sine kernel acting on a finite 
size interval [-q;q]. We determine its asymptotic resolvent as well as the first terms 
in the asymptotic expansion of its Fredholm determinant. Further, we apply our results 
to build the resolvent of truncated Wiener-Hopf operators generated by holomorphic 
symbols. Finally, the leading asymptotics of the Fredholm determinant allows us to 
establish the asymptotic estimates of certain oscillatory multidimensional coupled in- 
tegrals that appear in the study of correlation functions of quantum integrable models. 
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1 Introduction 



The sine kernel 

sin ^{A- u) 

S (A,ti) = — ^- (1.1) 

is a very important object in mathematical physics. In particular, the Fredholm determinant of 
the integral operator I - S acting on some interval J cR appears in random matrix theory [18]. 
In the bulk scaling Umit, dety [/ - 5] stands for the probability [19] that a matrix belonging to 
the Gaussian unitary ensemble has no eigenvalues in xJ. The kernel (1.1) also appears in the 
theory of quantum integrable systems. In particular, the determinant det/ [/ + y^], y being a 
parameter, describes various zero-temperature correlation functions of the impenetrable Bose 
gaz [36, 29]. 

In all these interpretations of the sine kernel, one is interested in the large x behavior of its 
Fredholm determinant. The first attempt to analyze the x +oo asymptotics of det/ [I - S] 
goes back to Gaudin and Mehta [18, 19]. In 1973, Des Cloizeaux and Mehta [15] showed that 

x2 1 

logdet[_i;i][/-5] = -y--logx + 0(l), x^+oo. (1.2) 

Three years later, using Widom's formula [40] for the asymptotics of Toeplitz determinants sup- 
ported on an arc, Dyson [16] gave a heuristic derivation of the constant terms cq and proposed a 
recursive method to compute the subleading coefficients ci, C2, . . . in the asymptotic expansion: 

logdet[_i;i][/-5] = -y-^logx + co + ^ + g + ... . (1.3) 

However, the forementioned results were somehow heuristic. It was only in 1994 that Widom 
[41] managed to prove rigorously the first term in the asymptotics (1.2): 

-ilogdet[_i;ij[/-5] = -^+o(l) . (1.4) 
dx 4 

One year later, this analysis was extended to the multiple interval case [42]. While Widom 
studied the asymptotic behavior of the Fredholm determinant by operator techniques, Deift, 
Its and Zhou applied the Riemann-Hilbert problem (RHP) formulation for integrable integral 
operators [27] to the sine kernel acting on a union of intervals UfJ^ and proved the existence of 
the asymptotic expansion (1.3). However, their method did not allow them to obtain an estimate 
for the constant cq, as they inferred the asymptotic expansion of log det/ [/ - 5] from that of 

P;, = x-^logdet[/-5] . (1.5) 
dx 

The first proofs of Dyson's heuristic formula for cq appeared in the independent, and based on 
two completely different methods, works of Ehrhardt [17] and Krasovsky [34]. 

We would like to point out that there is a very nice connection of the sine kernel to the 
Painleve V equation [29], as Px solves this equation. The link between Painleve V and Px 
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was also investigated in [13] in the framework of RHP. It was shown that one can deduce this 
Painleve equation directly from the RHP data. 



This article is devoted to the study a generalization of the sine kernel (1.1). This kernel, that 
we will refer to as the generalized sine kernel (GSK), is of the form 

V U,fi) = , / , [e^ (A) e- Ou) - e_ (A) e+ Ou) ], (1 .6) 

2i7T {A - n) 

where 

We will be more specific about the functions F, p and g later on. 

Various particular cases of the kernel (1.6) already appeared in the literature. These par- 
ticular kernels were mostly used for the description of correlation functions of matrix models 
or quantum integrable models equivalent to free fermions (see e.g. [39, 37, 24, 25, 32, 9, 10, 
26, 28, 8]). In the present paper we consider a rather general case, only based on the analytic 
properties of the functions F, p and g. The GSK (1.6) plays a crucial role in the study of cor- 
relation functions of (non free-fermionic) quantum integrable systems [31]. It is also useful for 
the asymptotic analysis of truncated Wiener-Hopf operators with Fischer-Hartwig singularities 
[33]. 

We investigate here the large x asymptotic behavior of the Fredholm determinant of the GSK 
in the framework of RHP. Our work is a natural extension of an unpublished analysis by Deift, 
Its and Zhou of the pure sine kernel / -i- by RHP. The pure sine kernel was also analyzed by 
RHP in [8]. 



This article is organized as follows. In Section 2, we announce the main results of the paper, 
namely, 

• the large x asymptotic behavior of the Fredholm determinant of the integral operator I+V, 
cf. (1.6); 

• the asymptotic resolvent of some Wiener-Hopf operators connected to (1.6); 

• the asymptotic behavior of coupled multiple integrals involving a cycle of kernels V (1.6) 
versus some holomorphic symmetric functions. 

The proof of the asymptotic behavior of log det[/ -i- V] is given in the core of the paper 
(Sections 3, 4, 5 and 6). More precisely, in Section 3, we recast the problem into a certain RHP. 
In Section 4, we transform this initial RHP into a RHP that can easily be solved asymptotically. 
This asymptotic solution is presented in Section 5 and used in Section 6 to obtain the leading 
and the first subleading terms of log det[/ -i- V] in the x — > -i-oo Umit. 

In Section 7, we apply these results to truncated Wiener-Hopf operators. We show how 
one can use the asymptotic resolvent of the generalized sine kernel to construct asymptotic 
resolvents of truncated Wiener-Hopf operators acting on [ -x ; x ], with x large. This asymptotic 
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resolvent is used to reproduce the low magnetic field behavior of the so-called dressed charge 
arising in the theory of quantum integrable models solvable by the Bethe ansatz [6]. 

Section 8 is devoted to the study of the asymptotic behavior of some particular type of 
coupled multiple integrals which can be obtained in terms of the GSK. This is in fact our main 
motivation to study the GSK: indeed, such integrals are closely related to those appearing in the 
description of correlation functions of quantum integrable models, and from the knowledge of 
the asymptotic behavior of the former one can obtain the asymptotic behavior of the latter, as it 
will be done in [31]. 

Finally, in Section 9, we consider the case of further modifications of the GSK, in particular 

those useful for the correlation functions of the integrable Heisenberg spin chains [31]. 

Some properties of confluent hypergeometric functions and proofs of several lemmas are 
gathered in the appendices. 

2 Problem to solve and main results 

2.1 Generalized sine kernel: assumptions and notations 

Let / + y be the integral operator with kernel (1.6) and acting on ])■ 

We assume that there exists some open relatively compact neighborhood U oi{-q;q\ such 
that the functions p, F and g, as well as the parameter y, satisfy the following properties: 

• F and g wee holomorphic on U, the closure of U ; 

• p is holomorphic and injective on U, p([-q;q]) c R, and p stabihzes the upper half 
plane "7/+ (resp. the lower half plane 'H-), i.e. piU n 9i±) c 9i±; 

• ye Do^r - [Ae C : \A\ < r}, where r is such that \rF\ < 1 and arg (1 + yF) e ] -;r ; ;r [ on 



We study the large x expansion of the Fredholm determinant of / + V under these assump- 
tions. This will be done by asymptotically solving a certain matrix RHP. It will become clear 
in the next section that the assumption p{[-q;q]) c R is tantamount to imposing the associ- 
ated RHP to be of oscillatory nature. Moreover, the case p{U C\ 'J-l±) c is obtained by the 
negation {y,g{A)) i-^ {-y,-g{X)). 

Before presenting the main result of this article, let us introduce some convenient notations. 
First, we define two auxiliary functions used all along the article: 



Note that /c is a function of the two parameters A and q, although we will sometimes omit the 
dependence on the second parameter. 



U. 



v{A)^-^\og{\+yF{A)) , 



(2.1) 




(2.2) 
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Finally, we will use the following simplified notations for the values of the functions p and 
y and of their derivatives at the points ±q: 

p±^p{X)\ , p'^=p'(A)\ , etc. (2.3) 

I A=±q I A=±q 

v±=vU)| , yl = v'U)| , etc. (2.4) 

I A=±q I A=±q 

2.2 The main results 

We now give the asymptotic behavior of the Fredholm determinant in the x — > +oo limit: 

Theorem 2.1. Let V be the GSK (1.6) with p, g, F and y satisfying the assumptions of Sec- 
tion 2.1. Then, in the x — > +oo limit, logdet[/ + V] behaves as 

log det[/ + y] = log det[/ + Vf^ + o(l), (2.5) 

with 

1 1 
log det[/ + V]^^^ = -ix J v{A)p'{A) dA - {vl + vl) \ogx- ^ v{A)g{A) dA 



+ log 



G(l,y+)G(l,y-)^^^(^;^) 1 f V{A)v{^D - v{A)V{p.) 

-2 -2 +0 d/ld/i , 

{2qp'J^{2qp'_y-Ky-{-q-q)\ 2J A-ix 



in which we have used the notations of Section 2.1. The Barnes G-function [3, 2] admits the 
integral representation: 

' . I.., ..A , .... r'(z) ^27) 



G(z+l) = (2;r)iexpJ-i^^^ + J?^/r(Od?L 5l(z) > -1, V^(z) = 





and we denote G(l,z) = G(l + z)G(l - z). 



r(z) 



Using the perturbation theory for singular integral equations one can refine the theorem and 

obtain sub-leading corrections. Although, in principle, nothing opposes to derive the next sub- 
leading corrections, the computations become more and more involved. We have proved the 
structure of the first corrections to the equation (2.6). 

Proposition 2.1. Let Vbe the GSK (1.6) with the conditions of Section 2.1. The leading asymp- 
totics logdet [/ -i- V]*^"^ o/logdet [/ -i- V] as defined in Theorem 2.1 has non-oscillating and os- 
cillating corrections. 

Let < 6 < q be such that the disks D±q^s of radius 6 centered at ±q fulfill D+q^s c U. Let 
s - 2sup^o^^uaD_g5 Then the first non-oscillating corrections are of the form 

^^Ol^i (2.8, 
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with 

The first oscillating corrections are of the form 

and the leading oscillating coefficient is given by 

Or = V y (^I'x^"^^^""-) ^i<rx(p.-P-) , (2.11) 



where we have introduced 

r(l + v+) \ K{q;q) 

«_ = e^^-"^ ^7^^^ {{2qp'_y-K{-q; q)}'^ . (2.13) 

r(i - V-) 

Remark 2.1. The GSK depends only on the combination ixp{A) + g{A) (see (1.7)). Therefore the 
Fredholm determinant and its asymptotics can only depend on this combination. This observa- 
tion allows us to obtain the complete asymptotic expansion depending on the function g(A) from 
the asymptotic expansion of the Fredholm determinant I + V corresponding to g = 0. Namely, it 
is enough to replace in the obtained formulae p{A) by p{A) - jgiA) and then expand into negative 
powers of x. 

It is quite interesting to apply the latter proposition in order to obtain the first few terms of 

the asymptotic expansion of det [/ + V]. The reason why we draw the reader's attention to these 
asymptotics is because they present a very interesting structure: the leading oscillating terms in 
the asymptotic expansion are just given by the sum of the leading asymptotics evaluated at v 
shifted by 1 or -1. This structure of the asymptotics seems to restore, at least partly, the original 
periodicity v — > v + n, n e Z, of the Fredholm determinant of / + y. 

Corollary 2.1. Let I+V be the GSK as above, det [/ + V]^'^^ [v] the leading asymptotics of its 

Fredholm determinant just as in Theorem 2.1, Ni and 0\ as in Proposition 2.1. Note that we 
have emphasized the structure of det [/ + [v] as a functional of v. Then the oscillating 
corrections 0\ can be reproduced from the non-oscillating part via the shift ofvby±\: 

det [/ + Vf^ [y] ^ = det [/ + vf\v+\] + det [I + Vf\v-\] . (2. 14) 

This structure of the first terms of the large x asymptotic expansion for det [/ + V] leads us 
to raise the following conjecture on the structure of the asymptotic series : 

Conjecture 2.1. The asymptotic expansion of the Fredholm determinant det [/ + V] of the GSK 
restores the periodicity v ^ v + n, n e Z, of the determinant. In particular, this asymp- 
totic expansion contains all the Z-periodized terms with respect to v of the leading asymptotics 
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det [/ + V]^^^ [v]. Thus, all the oscillating terms can be deduced from the non-oscillating ones. 
More precisely, let 

^ [V] (.) ~ det [7 + vr [V] (l + ^^"^ M + ■ ■ ■ ^ (^^f ) M ^ . . . ) (2.15) 

stand for the formal asymptotic series corresponding to the non-oscillating part of the asymp- 
totic series for log det [7 + V\ There Ck QQ [y] are polynomials of degree k in X whose coeffi- 
cients are functionals in v. Moreover each of the Ck 's has no oscillating exponents of the type 
Q±ixp±^ Yjjgfi the formal asymptotic series for det [7 + V] is given by 

det [7 + y] ~ ^ ^ [y + n] (x) . (2.16) 

n€Z 

This conjecture is supported by (2.14) and also by the results of [38] where several sub- 
leading corrections to the asymptotics of the Fredhohn determinant of the pure sine-kernel were 
computed. 

The first application of the asymptotic behavior of the GSK we consider in this article con- 
cerns the asymptotic inversion of truncated Wiener-Hopf operators. We will prove in Section 7 
the following proposition: 

Proposition 2.2. Let I+K be a truncated Wiener-Hopf operator on]-x;x[, acting on functions 
g e l2(R) as 

X 

{{I + K) .g] it)^git) + jKit-t')g (f) dt' . (2. 17) 

-X 

The kernel K is defined by its Fourier transform F: 

K(t)^r-'[F](t) , (2.18) 
and we suppose that there exists 6 > such that 

• F admits an analytic continuation to |z : |3(z)| < 5|; 

• ^ F ± i6) e (R); 

• the analytic continuation of\+F never vanishes for |3(z)| < <5. 
Then the resolvent 7-7? of I + K fulfills 



^ ^(^) |^e-(^-) - ^e-'>(^-)| . O (e--) , (2.19) 



logil+F(fx))^}. (2.20) 
p-A\ 




2inJ 
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Our main motivation to study the asymptotics of log det [/ + V] comes from the theory of 

one dimensional quantum integrable models. Indeed, the generating function of the zero tem- 
perature two-point correlation functions (at distance x) of different quantum integrable models 
[31] has a series expansion in terms of cycle integrals of the type 

In \rn\ y J ^2inf ^" {z} 1 ( _ ^^ ( _ ^ ^ ■ (2.21) 

There F ([ ; ^ ]) is a closed contour around [ ; ^ j whereas the function Tn is holomorphic in 
some open neighborhood of [ ; ^ and symmetric in the n variables {A} (we set An+\ = A\) 
and in the n variables 

In Section 8, using the above results for the GSK, we prove the following asymptotic ex- 
pansion of J,j Yfn^ in the x ^ +oo limit : 

Proposition 2.3. Let Tn and In [Tn] be as above. Then for x — > -t-cxj, 
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^" = 2^ {'^^'^^^ ^ ( [A . luV- 



e=0 



\og{2qp'^x))Tn 



{crqY 



(2^^^ J"" p{n-p){q-crA) 

n y I dAA^i i l{A + e], , {fif-P 

2 (2nf pj^in- p) {A - fi) f^^" \{A + e}, {Ar' , {/i}""^ 

_^^/{/^ + ^K{/^F-M^rMl (2 22) 



with 



(-1)"-! 5"!^ 

Cn - 



(-1)"-! a''logG(l,yo) 



(n-1)! 57" 

7=0 



(n-1)! dy" 

and where {A}" denotes the set formed by n copies of the same parameter A. 



, vo = — log(l+r), (2.23) 
7=0 27r 



Moreover, in Section 8 we will also describe the form of the sub-leading corrections to this 
result. 



2.3 Comparison with known results 

There are several results in the literature concerning the asymptotic behavior of the Fredholm 
determinant det [/ -l- yS]. This determinant corresponds to the GSK with p-iA,F-\ and g - 
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0. It is clear that we reproduce the answer concerning the leading asymptotics of det [/ + 75] 
analyzed in [7] and [4]. 

d 

As observed in [29], x — log det 1/ + I satisfies the fifth Painleve equation. The authors 
ax 

of [29] used this property to obtain an asymptotic expansion of log det [I + yS]. This fact was 
also exploited by the authors of [38] in order to derive the first few terms in the sub-leading 
asymptotics of the latter quantity. Their result reads 

d 

X— log det[/ + yS]- -4ixvo - 2^ - /— 
ax " X 



■2 ^'/r^'_„„^\2 /r^'_„„^\2 o-4!Jc 



+ ^ I (4x)^^o e^'-^ - — \ , (2.24) 



with vo given in (2.23) and ^ = 2. It is straightforward to see that in such a hmit A^^i = iv^ and 

0,^A_ U (2,.)-o iUl^]' . ^ (IMS] , (2.25) 
(2qf\ \Tivo)) (2^x)4^o \r(-yo)j J ' 

which reproduces the oscillating terms (2.24) after setting q - 2 and taking the q derivative. 



3 The initial Riemann-Hilbert problem 

The GSK (1.6) belongs to a special algebra of integral operators, the so-called integrable integral 
operators. This algebra was first singled out in [27] and then studied more thoroughly in [13]. 
It is well known that many properties of these integrable operators can be obtained from the 
solution of a certain RHP. 

In this section, we formulate our problem in terms of a RHP that we then asymptotically 
solve. 



3.1 Notations 

An important property of completely integrable integral operators is that their resolvent still lies 
in the same algebra. However, before presenting the formula for the resolvent we introduce 
some quite useful vector notations. Namely, let 

I (^) ) = (^)J , ^eHA)\= #(!)(-.- U) , (i)) , (3.1) 

so that the kernel V has a simple expression in terms of | (A) ) and ( (A) |: 

{E^{A)\E^Ui)} 

V(A,i^) = . (3.2) 

A- IX 

Observe that 

<£^(^)|£«(A)) = 0, (3.3) 
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and, hence, the kernel V is not singular at A- fi. 

Let I (A) ) be the solution to the integral equation: 

|F«0")) + J yU,/x)|F«U))d^ = |£^0")), (3.4) 

and ( (A) \ be the solution to the corresponding dual equation. It is convenient to write 
I (A) ) as well as its dual < (A) | in a form similar to | £^ (A) ) and < £^ (A) |: 

I /r« (^) ) = yjm. |A (^)J , < /r^ (^) I . VFU) (-/- (A) , U U)) . (3.5) 

Then the resolvent reads: 

, {F^{A)\F''{^)) yJE^E^..,,..,. 

R{A,n) = = — U+{A)f-{ix) - f+(j/)f-(A)\ . (3.6) 

A- 1^ 2m {A - n) 

3.2 The Riemann-Hilbert problem associated to a GSK 

Proposition 3.1. Let V be the GSK (1.6) understood as acting on ([ -q ;q]), and such that 
det [/ + y] ^ 0. Then, there exists a 2 X 2 matrix x {A) such that 

\F''{A))^X iA) I {A) ), < (^) I = < F^ (^) Uf-i (^) . (3.7) 
The matrix X {A) is the unique solution of the RHP: 

• X is analytic onC\[-q;q\ ; 

• X+(X)Gx{A)^X-iX) forAe\-q;q[ . 
The jump matrix G^for this RHP reads 

Finally, x and its inverse can be expressed in terms of\F^ (A) ) and of its dual { F^ (A) \: 

, r\F^M)XE^ _i , , f |FV))(FV)| ^ 
^(i) = /2 - dju, ^ (i) = /2 + dju. (3.9) 
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Figure 1: Original contour for the RHP. 



We emphasize that the big O symbol, ^ | | | j, is to be understood entrywise. Moreover, 

X± (m) stands for the non-tangential limit ofxi^) when A approaches a point /x belonging to the 
jump curve from the left, resp. right, side of the contour (see Fig. 1). 

Proof — The unicity of solutions to this RHP is proved along the same line as in [35]. The 
proof of existence of solutions is based on the equivalence between RHP and singular integral 
equations which, in the case of the above RHP, implies 

1 

XU)^I2 + J ^X+(M)\E''(M)){E''(4i)\ ,^eC\[-^;^J. (3.10) 
-g 

The solution to this equation can be expressed in terms of the resolvent kernel I - Rof I + V 

1 

X(A)^l2 + J ■^{\E'').(I-R)}(M){EHfi)\. (3.11) 
-g 

In its turn, the resolvent kernel exists as det [I + V] 0. Moreover, the explicit construction 
of the resolvent through a Fredholm series shows that (/i,/i) ^ /?(/l,//) is analytic in U x U. 
Hence, so is |F«Ou)) = [\E'^).(I - R)\(ju). The estimate \x\ = o(log _ ^2|J ^ ^ ^ ±q 
follows from the integral representation (3.11) supplemented with the fact that both ( \ and 
I F'^ ) are smooth on [-q',q]- 

Applying (3.9) to lE'^iA)} and {E^{A)\ we obtain the equations (3.7). Hereby one can easily 
check that due to the orthogonality condition (3.3) the transform (3.7) is continuous across 
[-q,q]. □ 

It is also possible to express logarithmic derivatives of det [/ + V] either in terms of the 
resolvent /? of / + y or in terms of;^'. Indeed, we have the 

Lemma 3.1. The derivative of log det [/ + V] with respect to x is related to the following trace 
involving the matrix x 

5;,logdet[/ + V] - d) —p{A)\i[d;oci^)o-3X~'i^)], (3.12) 
r([ -?;?]) 

with (T-i = ^\ whereas its derivatives with respect to y and q are expressed in terms of the 
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resolvent as 

1 

5ylogdet[/ + y] = J—R(A,A), dglog det [I + V] ^ R(q,q) + R(-q,-q) . (3.13) 
-g 

Proof — The last two equations are easily proved by the multiple integral series expansion of 
log det [/ + V]. We shall only focus on the equation relating the x derivative of log det [/ + V] 
tox- Clearly 



log det [/ + V] - J [d,V. (7 - R)] {A, A) dA , 

-9 



with 



d^V(A. 



r([-9;?]) 



(3.14) 



(3.15) 



So that, using the representation (3.6) of the resolvent R in terms of ( | and | ) and the fact 
that < F^iii) I era I F\A) ) - tr [(t^\ F\A) ){ F^(n) |], we get 



..r.^T/i ^ f ^^ (£^(i) |cr3|£^(i)) 
5.1ogdet[/+y] = - -piz)j dA 

r([ -?;?]) -9 



+ tr-; 



^^/^(z) J did/i|F«(i))<£^(^) 



1 1 \ |£«0u))(F^0u) 
X 1^ T 0-3 



A- z A- ji 

Using the integral expressions (3.9) for x and;^f~\ we obtain 

{E^{A) \o-^\E^{A)) 

- py^Z) I u/i 

r([-^;?]) 



(3.16) 



5;,logdet[/ + V] = - ^^/'(z) J di- 



+ tr 



lr([-^;^]) 



= ^ ^p(z)tT{d^(z)(r3X-'(z)}, 



(3.17) 



r([ -«;«]) 



where we used (3.7). 
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It is worth noticing that formula (3.12) is particularly effective when p is a rational function 
as then the contour of integration can be deformed to the poles of p (including the pole at oo). 
The integrals can be then easily calculated. In particular, in the case p {A) = A, we have the 
following result: 

Corollary 3.1. Let Xi be the first non-trivial coefficient of the expansion of x around oo, i.e. 
x{X)^h + ^+o(-\. (3.18) 



A \A^ 
Then 

log det [/ + V] |p=id= -^tr {xio-3] . (3.19) 

Proof tr (0-3;^^ 1 ) is the residue of the pole at infinity of /I i-> /I tr {^^(/l) o-^x' ^W]- ^ 

In this way, we recover one of the formulae derived for the pure sine kernel [13], but also 
for more general kernels as in [24, 25, 28]. We emphasize that (3.19) is valid not only for the 
pure sine kernel as it was originally derived, but also for the generalized sine kernel with p = id. 



4 Transformations of the original RHP 

In this section we perform several transformations on the RHP for x so as to implement Deift- 
Zhou's steepest descent method [14]. The first substitution maps the RHP for the matrix into 
a RHP for a matrix 3 whose jump matrix has 1 on its lower diagonal entry. This jump matrix is 
then easily factorized into upper/lower triangular matrices. This factorization allows us to define 
another RHP for an unknown matrix T whose jump matrices are already exponentially close to 
identity uniformly away from the endpoints ±q. It remains to construct the parametrices at q 
and -q. These parametrices enable us to define a matrix n satisfying a RHP with jump matrices 
uniformly I2 + (1) when x — > +00. 



4.1 The first step x ^ 

Let 



a (A) - exp 



/ 



p. - A 



dju 



^K{A) 



A - q 
A + q 



viX) 



Then clearly, a {A) solves the scalar RHP 

a-{X) - a+{A){\ +yF {A)), Ae[-q;q\, q;(/1) — > 1 at /I — > 00 



(4.1) 



(4.2) 



The functions k{A) and v{X) were already introduced in (2.2) and (2.1). In the following, we 
shall also use another representation for the function a{A) : 



a {A) = Kp{A) 



piA) - piq) 



p{A) - p{-q) 



v{X) 



(4.3) 
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where Kp is defined as 



\ogKp{X;q) = \ogKp{X)^ ffvU) -^)d/x. (4.4) 

J \ p{A)-p(jj) A-jjL) 
-q 

We specify that we chose the principal branch of the logarithm, i.e. arg e ] -;r ; ;r [. Due to our 
assumptions on y, F and p, Morera's theorem implies that the functions v, log/f and \og Kp are 
holomorphic on U . Moreover we have |!R(y(/l))| < 1/2, V/l e U. Indeed 

y (i) = log 1 1 + yF (i)| - J-arg (1 + yF {A)) , (4.5) 

and we have assumed that arg (1 + yF) e ] -;r ; ;r [. 

We use the function a to transform the RHP ioxx- Let us define the matrix 5 (/I) according 

to 



sw)=xw)r« (4.6) 



This new matrix H {X) satisfies the following RHP: 
• 3 is analytic on C \ [ ; ^ ]; 



. |H(^)| = O ( I I j _ ^2p-3^R(v.) _ for ^ ^ ±^ ; 

. S(^) ^ /2 ; 

. H+(^) GH(i) = 5_(i) for i e [ ; ^ ] . 

Here the new jump matrix Gh reads 

/1+P(i)e(i) P(^)e«P(^) \ 
"^ ^"i fi(i)e--''(^> 1 j' ^^-^^ 

and 

= TT^^'W''"' = ^e«"'. (4.8) 

1 + yF{?C) a\ {A) 

= --^-^^ - 2.-e-(^) ^^^^ a! (i) . (4.9) 

1 + yF{A) qsW 

The solution of this RHP for S exists as it can be constructed from^. Moreover it is unique 
as seen by arguments similar to those providing uniqueness of the solution to the RHP for;^^. 
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4.2 The second step S ^ Y 



As already mentioned, the jump matiix Gn admits an explicit factorization into a product of 
upper and lower triangular matrices: 

Gh = M+M_. (4.10) 

The matrices M± are given by 

/ 1 P ( A^ e'^P'-'^^ \ / 1 \ 

M.U) = (o '\ ], M_U).(^^^^^_,,,,, J, (4.11) 

and can be continued to J7 n ']-(+, resp. U n H-, where we recall that 'J-(± is the upper/lower half 
plane and U is the domain of holomorphy of all the functions appearing in the RHP. Then we 
draw two new contours r± in ;? (?7) and define a new matrix T(A) according to Fig. 2. 




Figure 2: Contours r+ and r_ associated with the RHP for Y. 



As readily checked, Y (A) is continuous across ] ; ^ [ and thus holomorphic in the interior 
of r+ U r . We have thus removed the cut along [ ; ^ ] and replaced it with cuts along r+ U r_ . 
The matrix Y solves the following RHP: 

• Y is analytic in C \ r+ U r_ ; 

• Y(^) - O 



Y(^) = O 
Y(i) = O 





\A + 







\A + 













\A + 





A — > ±q; 



l)|i + #'^^^(^-Hog|i + ^|, A -j> ±q; 

1 / AeHjii 
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A—>oo 



Y+ (i)M+(i) = T_ (A) 
Y+ (i)MziU) - T_ (A) 



for /I e r+ , 
for ^ e r_ , 



where the domains Hj, Hu, Hm are shown on the Figure 2. 

Clearly, the solution of the RHP for Y exists and is unique. Hence, the matrices Y andx are 
in a one to one correspondence. 

Note that, except in some vicinities of q and -q, the jump matrices M+ and MZ^ for Y are 
exponentially close to the identity matrix. Therefore, to study the asymptotic solution of the 
RHP, it is enough to study the local problems in the vicinities of q and -q. 



4.3 Parametrix around -q 

We first present the parametrix on a small disk £>-^,<j c ?7 of radius 5 and centered at -q, that 
is an exact solution of the RHP: 



f is analytic on £>-^,<5 \ {r+ U r_} ; 
PiA) = O 



P{A) = O 
P{A) = O 



+ 



log \A + q\, A—^ -q; 



log \A + q\, A -q; 

X&Hu 



\\\A + ^r'^32l(v-) log 1^ + ^1 ^ A ^ -q; 



uniformly for A e dD-q^s, 



y\—e 



P+(A)M^iA) = P^(A) for i e r+ n D_^,^, 
P+iA)MzH^) ^P-U) for ^ e r_ n D_^,5. 

Here s = 2sup|5l(v)| < 1. The canonically oriented contour dD-q^ is depicted in Fig. 3. 

The RHP for f admits a class of solutions. Each element of this class is related to another 
one through a left multiplication by a holomorphic matrix that is uniformly h + o{\ lx^~'^^ on 
dD-q^s- In order to construct the solution P to this problem, we first focus on the simpler case 
where the functions F, g and Kp are constant. Then the solution to the RHP for Pconst can be 
obtained by the differential equation method [23, 11, 12]. This leads to the solution 



Pconst (/I) ^'¥(A)L (A) [-/^-J"'"' . 
Here = x{p{A)-p {-q)), y = / log (1 + jF) jln. 



»P (- V, 1 ; -i^-q) /Z712 ^ ( 1 + y, 1 ; K-q)\ 
-ib2i ^ (l - y, 1; -i^-q) ^ (y, 1; i^-q) J 



(4.12) 



(4.13) 
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Figure 3: Contours in the RHP for P. 



and finally 



sin[;rv]r-(l + v) 

2v' 



Jxp(-q)+g 



7rK^[x(p(q)-p(A))] 



nKi[x(p(q)-p(A))f'' . . , 
b2i (A) = -i ^ . , e-'^Pi-D-g 



(4.14) 



(4.15) 



sin [nv] (v) 

^(a,c; z) denotes Tricomi confluent hypergeometric function (CHF) of the second kind (see 
Appendix A). It solves the differential equation 

zy" + (c-z)y'-ay^O. (4.16) 

Recall that ^ has a cut along R.~. Note that this choice for the cut of implies the use of the 
principal branch of the logarithm: -n < arg (z) < n. The expression for the piecewise constant 
matrix L depends on the region of the complex plane. Namely, 



L(A)^ 





g-2OTy 
Uttv q 







1 



-n/2 < arg [p (A) - p {-q)] < n/2, 
nil < arg [p {A) - p {-q)] < n, 

-n < arg [p (A) - p (-q)] < -n/2. 



(4.17) 



The reader can check using the monodromy properties of Tricomi CHF (A.4) and (A.5) that 
the jump condition for constant functions F and g are satisfied by the matrix P const- Moreover 
the asymptotic expansion for *P (a, c; z) at z ^ oo allows one to check that Pconst has the correct 
behavior at infinity. 

In order to extend this result to the case of arbitrary holomorphic functions F (A), g (A) and 
Kp (A), it is enough to add the A dependency in all places where these functions appear. One 
ends up with the following solution to the RHP for P : 



P(A)^^(A)L(A)[-i^-,] 



e 2 



(4.18) 
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Here ^x(p(A)-p (-q)), 



^{-y{A),\--i^^,) ibu(A)^{i + y{A),l;i^-,} 
-ib2i (i) »P ( 1 - y (i) , 1 ; -K-q) *F (y (i) , 1 ; J 



(4.19) 



with 



• sin[;ry(/l)]r2(l +y(^)) „(X)+ixv(-q) )^ ^ 

bn U) = ; T^TTTTS^ ^ ^ = -IV {A) u {A; x) , 



7T4(A)[x(p(q)-p(A))] 



2v(A) 



(4.20) 



sm[;ry(/l)]r^(y(/l)) u{A;x) 

and finally 

« ^) - rn^'S! {'^^('^) ^^^'^ ^^^^-P ^^^^^'V (4.22) 

In the above formulae we have explicitly stressed the dependence of the functions bn, t>2i 
and y on A. Finally, the matrix L (A) is given by (4. 17) with y replaced by the function y (A). 

This construction originates from the observation that the replacements F F {A), g i-> 
g (A) and Kp Kp (A) preserve the jump conditions as the latter hold pointwise. Of course, 
once the parametrix P is guessed it is not a problem to check directly that it solves the RHP 
in question. The asymptotic behavior is inferred from (A.6), whereas the jump conditions can 
be verified thanks to (A.4) and (A.5). Furthermore, due to the definition of the matrix L, the 
solution is continuous across the line arg [p{A) - p {-q)] = n and thus analytic in the whole 
domain {A e C;%[p(A) - p (-q)] < 0}. 

4.4 Parametrix around q 

The RHP for the parametrix P around q reads 

• ^ is analytic on ^ \ {r+ U r_} ; 

• P{A)^ol] ]]\A-qr^^'-hog\A-q\, A —> q ; 

• p(A) = 72 + uniformly for A e dDg^g; 

^ r ^+ (A) M+ {A) = P- {A) for ^ e r+ n D^,^, 

• I ^+ {A) Mil {A) = P. {A) for ^ e r_ n D^^s. 
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Figure 4: Contours in the RHP for P. 



and e = 2 sup^^,^^ < 1 . The solution to the RHP for the parametrix P around q can be 

formally obtained from the one at -q through the transformation q ^ -q and v — > -y on the 
solution to the RHP for P. Indeed, the two RHP are identical modulo this negation. 
Just as for the parametrix around -q, we focus on the solution 



where ^q = x{p {A) - p (q)], and 

w{v{A),l- -i^,) ibu (i) ^ (l - y (/I) , 1 ; i^,) 

-ib2i {A) >F (l + y U) , 1; -i^q) »P (-y U) , 1; /^J 



(4.23) 



(4.24) 



Here 



hniA) . , sinkv(^)] T\l - y(^)) ^^^^^^ _ ^ ^^^^^^ 

TTK^iA) 

l>2l(A) - i :^ TTT7 - ' , 

sin[;ry(i)] r\-v{A)) [x{p{A) - p(-^))]2>'(-*) u{A; x) 



and 



u(A;x) - 



r(l-y(i)) \ x'^^>[p{A)-p{-q)] 

r(i + y(^)) 



v{A) 



Jxp++g{A) 



KpiA) 



(4.25) 



Just as for the parametrix around -q, the matrix L (A) depends on the quadrant of the complex 
plane: 

^ h -nil < arg [p {A) - p {q)] < njl. 



L(A) = 



( I qiZu) ) <arg[p(A)-p (q)] < n, 

-n < arg [p (A) - p (q)] < -n/l. 



(4.26) 



g2OTv(/i) Q 

1 
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4.5 The last transformation T ^ 11 



Let 

U(A)^\ r(A)P-\A) for A eD^q^s, (4.27) 
TiA) forAeC\[Dg^s^D_g^s}. 

Introduce the curve = {r+ U r_)n|D^^5 U £>-^,<j}. Then His continuous across '^\{^, -q}. 
Since 11 is holomorphic in a vicinity of ^, we have that 11 is holomorphic in Dq^g U D-q^g \ 
[q, -q]. This, in turn, due to the estimates for P, P and Y around the points ±q, ensures that 
the singularities at these points are of a removable type. Hence 11 is holomorphic on the disks 

(5 U D-q^s. Finally, we see that 11 satisfies the following RHP: 

• n is analytic in C \ En (c/ Fig. 5) ; 

• nU) = h + 0{\IA) for A^oo; 

{ n+ {A) M+ (A) - n_ (A) for AeT'^, 

n+ (i) MI = n_ (A) for Aer_, 

n+ (A) P (A) = n_ (A) for A e dD.q^s, 

\ n+ {A) V {A) = n_ {A) for A e dDq^g. 

The solution to the RHP for H, exits and is unique as seen by standard arguments. 







En = r: u r; u dD_q,s u dOq^s 



Figure 5: Contour En appearing in the RHP for H. 



The jump matrices for H are uniformly exponentially close to I2 in x on Fl U F^ and uni- 
formly I2 + 0(x^^') on dDq^s U dD-q,6^ with s - 2sup5£,^^uaD_,i As a consequence, h 
is the unique solution of the RHP, up to uniformly O (^^^') corrections. In addition, using the 
equivalence between singular integral equations and RHP, the asymptotic expansion of H can 
be obtained by a Neumann series. This will be done in the upcoming section. 
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5 Asymptotic solution of the RHP 

In this Section we asymptotically solve the above RHP for n. 

We derive an asymptotic expansion into negative powers of x for the jump matrices for H, 
and use it to prove the existence of an asymptotic series for IT. The corresponding asymptotic 
series for;^' follows readily. One can finally infer the asymptotic behavior of the resolvent of the 
GSK up to any order in 1 jx. 

5.1 Asymptotics of the jump matrices 

Denote the jump matrices for 11 by /2+A (/I). Then the matrix A {X) has the asymptotic expansion 
in the Umit x — > +00: 

m=l 

with e = 2^w^aD,^^dD.,^ 

The explicit form of the matrices A*^") (i; x) depends on the position of A in the contour 
En: they vanish to any order in l/;c on r+ U Fl, whereas the asymptotic expansion for A on 
dDq^s U ^-q,s follows promptly from the asymptotic expansion of Tricomi CHF (A.6). More 
explicitly, for any n e N*, 



A(") {A; x) = 



for /I e dD-q^s, 



A[_tU;x) 
{p{A)-p{-q)\ 

a["^(A;x) 

for A e dDn §, 

[p{A)-p{q)f '''' 

for^er;ur 



(5.2) 



We have separated the jump matrices into their pole parts [p(A)- p±] " and regular parts a["^ 

n u(A; x) \ 



and a[J, with 



A^]iA;x)='- 
* ' n! 

for A e dD-q s, and 

aS(^;x) = ^ 



1 

n 



v{A)u{A; x) 



v{A)u{A; x) 



viA) 
1 



n u{A; x) \ 

v(^) 
1 



{vml 



(5.3) 



{-\f{v{A))l 







(5.4) 



for A e dDqs . Here we use the standard notation (y)„ = r(v + n)/r(y) and uiA;x), resp. 
u (A; x), have been defined in (4.22), resp. (4.25). Thus, the matrices A*^"^ sUghtly depend on x, 
but their entries are a 0(;c^). 
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5.2 Asymptotic expansion for 11 

Using the equivalence between RHP and singular integral equations we can express IT in terms 
of its boundary value from the "+" side of the contour Xn 



n(/l) = /2 + 



2n 



d^ 
A - s 



(5.5) 



In its turn Ti+(A) belongs to (Sn) and fulfills the linear singular integral equation of Cauchy 
type 



n+(z) = /2 + cj„[n+A](z). 

Recall that the Cauchy operator on (En) is defined as 



(5.6) 



^ C2 llgll, where ||.|| is the canonical L?- (En) 



C^nfeKz) = l_im Csn[g](0 and C^^lgM = — J ^IJ—, t^Tn. (5.7) 

The notation t ^ z'^ stands for the non-tangential limit of t approaching z from the " + " side of 
the contour En. Recall that the Cauchy operator is bounded: i.e. there exists a constant C2 such 
that, for any function gel? (En), one has CJ^ [g] 
norm. 

The matrix 11+ can be asymptotically approximated by the following series: 
Proposition 5.1. Let 11+^ be defined recursively according to 

nf = J]c^^[nf^^A(^)] with nf = /2. 

p=i 

Then, for any integer M > 0, there exists a constant C(M) > such that 

M-\ 



(5.8) 



p=0 



C(M) 



(5.9) 



Proof — Let us prove this statement by induction on M. For M = 1 we have that 

lin+ - /2II = ||c+^ [(n+ - h) A] + cl^ [A] II < C2 ||n+ - /2II i|A|| + C2 ||A|| . 

Therefore, for x large. 



Iin+-/2||< 



C2IIAII C(l) 



I-C2IIAII - X 



l-e ■ 



(5.10) 



(5.11) 
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Let us now suppose that the result holds up to M. Then, 



M 






M 




k=0 






k=l 











M 


( M- 


P ^ 
























p=i 


k k=i 


) 





M 



C(M + 1) 



(5.12) 



for some constants Ca^, and C(M + 1). We used the fact that all A^^ are in (Sn) and that 
||n+|| is bounded in virtue of (5.11). □ 

Let us now extend this result for points A being uniformly away from the contour Sn. Define 
the matrices 



n(0)(z) = l2, U^P\Z) = 2 K-'^A(^)] (Z), > 0, 



M 



n(z;M) = J]x-pn(p>(z), 

p=0 

which are analytic away from Sn. Then we have the following result: 
Proposition 5.2. Let K be any compact subset of C\ Sn. Then, VA: e N, VM e N*, 

kl C{M) Igth(Sn) 



(5.13) 
(5.14) 



\d\ll{X)-d\ll{A;M-\)\ 



AeK. 



(5.15) 



i^ere |.| denotes the usual max norm |n| = max; y d(., .) is any distance on C and Igth(En) 
is the length of the curve Xn. 



Proof — Let /t e N, M e N*, then 

M-l 



d\U{A)-Y,^~'d\Yl^''\^) 



p=0 



2inJ 



ds 



M-l 



(A-s) 



k+l 



I n^(^) A(^) - 2 Z "^""''(^^ I 



/c!C(M)lgth(Xn) 



due to (5.9). 
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5.3 The functions /± to the leading order 

We now perform the transformations from n back tox- 
The solution to the RHP of Proposition 3.1 reads 

XiA) = n{A)x'''\A). (5.16) 

We call;^f(*'^ the zero* order solution (i.e. obtained for 11 = I2). In the vicinities of the endpoints 

of [ ; q],x'"^^ is given as 



/'\A) = 



M;1 (A) aiA)--^^ , A e n {0 < sxg[p(A) - p-] < n/2} , 
PiA) M;1 (A) a(A)-^^ , A e D^^s n {nil < axg[p{A) - p+] < n] . 



Similarly, on [ ; and uniformly away from the endpoints, 

X^'^KA) = M-\A) a^iA)-"' , Ae\-q + S;q-6[ . 
In the 3 (A) = 0+ Umit and for %iA) e [-q;q], 

so that, for Ae]-q + 6;q + 6[, 

( (A) ) _ jTTviX) r„.-l /- 1\ „ ^^ j7n'(X)Vi ( 1 



fPiA) 



j = e'-('^)[a-i(^)e,(^)e'-(-^)f ( J ), (5.18) 



where we have explicitly written all the dependence on A. 

When A e [ ; + <5], we should multiply the latter expression by P. Using the decom- 
position (A.7) of Humbert CHF into a sum of two Tricomi CHF we get 

( fr^iA)]_^.^( e^A) r(l + v)0(-v,l;-/^_,) \ 

iyl°^U)J Up(A)d [ r(l-v)cD(v,l;/^_,) J' ^^-'^^ 

with ^q = x[p (q) - p (A)] and = x[p (A) - p (-q)]. 
Analogously, for Ae [q- 6;q], 

(fi'hA)] .^(e^\-U r(l-v)0(v,l;/^,) ] 

[f^'\A))-' [ K,iA) ] [r(i+v)a>(-v,i;-.-^,) J- ^'-'"^ 

Note that the piecewise expressions for the functions f^^ {A) are in fact analytic in a vicinity 
of their respective domain of validity, although they have been obtained by taking the Umit of A 
approaching a point in [ ; ^] from the upper half plane. More precisely, the formula (5.19) 
holds on D-q^s, (5.20) on Dq^g, and (5.18) on the connected component of the interior of !,„ 
containing [6 - q;q - S]. This observation follows from (3.3), but of course it can be checked 
by a direct computation based on the expression for the matrix in the lower half plane. 
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5.4 Integral bounds for the resolvent. 

We now introduce a function R^^^ and show that it is a good approximation of the resolvent 
in the sense that 

tr(i?-/?(°^) = 0(x^-i) . (5.21) 
Such estimates are necessary for the integration of the y-derivative of log det [/ + V]. 

Definition 5.1. Let 11^ (A) denote the solution of the RHP given in Subsection 4.5 whose jumps 

are on circles of radius t and on the corresponding curves Y'^ and Y'_. 

We can then write the solution of the RHP for;^' as;^f (i) = HtX^^ ■ Yhere x^? do not depend 
explicitly on t. The radius t only determines which patch we should use for the definition of 
the matrix Moreover the whole combination YIjX^^^ does not depend on the radius r at all. 
Hence, we can represent the exact resolvent as 

R iA,p) ^ ( (A) I [x?^ (A)]-' ^^IS^^lMH^m (^) I (^) ) . (5.22) 
L J A- p 

There, without altering the value of R (A,p), we can chose different values of t depending on 
the point we consider. This is quite useful as we can take one value of the radius r in order to 
have estimates around ±q and another one to perform estimates in the bulk [6 - q;q - 6]. This 
will become clearer during the proof of the proposition below. 

Definition 5.2. Let us fix 6, q > 6 > and define what we call the diagonal zero^^ order 
resolvent 

R''' U, A) = ^ (d,ff (A) ff^ (A) - d,f!_'^ (A) (A)) , (5.23) 

where the functions f^\A) are given by (5.18) for A e [6 - q;q - 6], (5.19) for A e [-q',S - q[ 
and (5.20) for A€]q - 6;q]. Similarly, \ F^'(^) (A) ) and { F^''^^ (A) \ are defined in terms of the 

r ■ AO) 

same functions f^ . 

We stress that the radius r previously introduced to build the exact solution IIt{A)x'^^ (A) 
and 6 appearing in the definition are, a priori, unrelated. 

Proposition 5.3. Let R(A,p) be the exact resolvent of the generalized sine kernel. Then 

tr(/?-/?(°)) = 0(x^-'), (5.24) 
where the O is uniform in y e Z)o,r- 

Proof — According to the preceding observations we have, for /I e [-q;-q + 6[u]q - 6;q], 
R (A, A) = R^^^ (A, A) + { F^-("> (A) \U2s (A) d;iU2, (A) \ F^'^^^ (A) ), (5.25) 

and 

R (A, A) = (A, A) + < (A) \Us/2 (A) ^;^Us/2 (A) \ F«'^«) (A) ), (5.26) 
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for A e [6 - q;q - 6]. The advantage of using two different matrices 11 for the corrections 
of R (A, A) with respect to the zero* order resolvent (A, A) is that the corrections are always 
analytic on the whole domain where they are considered. One does not need to take into account 
that (A) has a jump across A - ±iq - 6). This might be problematic as, for instance, the 
integral of 5^11 (A) on [-q;6 -q] might be ill-defined. Moreover the uniform estimates that 
we have derived for the matrix IT (A) only hold uniformly away from the jump contour. As we 
will only integrate the terms containing 1125 on [-q;6 - q[u ]q - 6;q], we will be in this 
situation. The same holds for the terms involving Ils/2- However, we would not be able to 
use the uniform estimates (5.2) for S^IT^ when integrating it on [ ; 5 - ^], as we would not 
always be uniformly away from the boundary of the jump contour for 11$. 
With this way of understanding the corrections we have 



6-q q 

tr (r - = U + J jd^ < f (A) \U2s (A) d^U26 (A) \ F'''^''^ (A) ] 

q-6 

J dA { (A) \Us/2 (A) dxUs/2 (A) \ (A) ). (5.27) 



6-q q 

\ .... 

q-6 

q-6 

+ 
6-q 

Let us start by the bulk part of integral, i.e. the part on [6 - q;q - 6]. From the expUcit 

form for f^^ on[6 - q ;q - S] given in (5.18) we see that these functions are uniformly O (1). 
Moreover, the uniform estimates for the matrices 115/2 (A) for A uniformly away from the jump 
contour guarantee that 

< F^'(O) (A) \Us/2 (A) d;fls/2 (A) I F^'^""' (A))^0 (/"i) , (5.28) 

the O (■^^^) being uniform in y, at least for y small enough. 

The situation at the boundaries is a little more complex. We only consider the right boundary 
as the other case is treated similarly. We still have that 1125 (A) = I2 + O [x"^'^) and 5^1125 (A) - 

O (.^'^ ') uniformly on [q- S;q]. However the functions f^^ (A) are no longer uniformly a 
O (1) on this interval. We should thus estimate the following integral 

1 

f f^^\^)ff(^)G^,^'(A)dA (5.29) 

with Gcr,cr' (A) = O being related to the entries of 115/2 (A) dAns/2 (A). The situation being 
similar for all the possible choices of cr and cr', we explain the mechanism for (cr, cK) = (+, +). 
The asymptotics of Humbert CHF guarantees that 

a)(a,l;±«Y) = ^(l + o(l)) (5.30) 

for some computable constants c± depending on a. These constants are continuous with respect 
to a belonging to an open neighborhood of v ([ ^ - 5 ; ^ ]), and so is the o (1) term. Hence, there 
exist an a independent constant C such that 

(l + |?|)^(«)a)(a,l;±jO <C . (5.31) 



28 



Indeed the latter function is continuous on R and has a finite Umit at oo. Moreover the constant 

C can be chosen in such a way that the estimate holds for a belonging to some small vicinity of 
v([q - 6;q]). Hence, by exphcitly extracting the x^~^ factor coming from G+,+ (A) we get that, 
for some constant C, 

<C'x^-'<p,ip{A)-piq)) , (5.32) 

with (f^ (t) = ^^^(''('i)) (1 + X The function cp^ (t) fulfills 

\ip, ip (A) - p iq))\ <C\piA)-p (^)r22^W^» (5.33) 
as, for any a e R, f i-^ t" / {\ + t") h bounded. The latter function is integrable on [q - S;q] 




One can prove, in a very similar way, the estimates for the Hilbert-Schmidt norm of the 
resolvent. Namely, 

Proposition 5.4. Under the assumptions of the previous proposition, 

\\R-R^%^0(x^-^) (5.34) 
with II.II2 being the Hilbert-Schmidt norm. 



/ru)/ru)G+.+(i) 



5.5 Asymptotic expansion of the resolvent 

We now prove that the asymptotic expansion for 11 can be used to obtain an asymptotics expan- 
sion for the diagonal of the resolvent R(A, A). We derive point-wise bounds for the latter as this 
quantity appears in the ^-derivative of the Fredholm determinant: 

dq log det [/ + y] = /? (q, q) + R (-q, -q) . (5.35) 

We need to estimate the error when we replace the exact resolvent R by the approximate one 
The magnitude for the error term follows from the following result: 

Proposition 5.5. Letx^^^ be the solution of the RHP for x up to the leading order in x, that is to 
say the one obtained from 11 = /2 and corresponding to the contour En with disks D^^ g having 
radius 5. Define the leading vectors { F^'^^^ \ and \ F^'^^^ ) as 

{ (A)\^{ £^ (A) U^°>(i)-i , I (A) ) = Ar^o>(^) I (A) ) , (5.36) 

and the leading order of the resolvent by 

«'0>(-^.).<^""'<-'>l^"°'»». (5.37) 

A- p 



Then 



^ R^pHA,A) 
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for X uniformly away from Sn and belonging to {-q;q \ Here, 

R^^'> (A, A) = -<F^;«» (A) I (A)}, (5.39) 

R^^ (A, A) = -{ (A) I U^p\A) I F'^'^^^ (A) ), p > 0, (5.40) 

in which 

(^; /:) d^n (A; /:) = g n^^^i) x"^ + O {^■^jjj^'^ ■ (5.41) 
Proof — Clearly, 

A- ju A- n 

^ -{F^'(0\A) I d;iF'^''^^\A)} - { F^'^^\A) \n-\A) dM^)\ Z^^'^^V) )• 

/l-»jU 

The corrections to the leading order for the resolvent are given here by the second term. 

The inversion operator on M2 (C): u is continuously difFerentiable around the iden- 

tity h- Thus there exists an open neighborhood W of the identity matrix I2 and a constant C > 
such that, V A, B e W, one has - < C ||A - Here ||.|| denotes any matrix norm. 
The matrices IT (A) and n {A; k) belong to W for x sufficiently large, as they both go to h in the 
X — > +00 limit for A uniformly away from En, and we get, from Proposition 5.2, 

||n"^ {A) dxTV {A) - IT^ {A; k) Oa^ {A\ k)\ 

<c||n(i)-n(^;^)llll5^n(^)ll 

+ c (iin (i)|| + iin (i) - n {a- k)\\) {A) - {A- k)\\ 

^ C{k) 

for some constant C{k). Thus, uniformly away from Sn and on the real axis, one has 

|(F«0)(i)|n-i(i)5an(^)-n-i(i;^)c^4n(i;/:)|F«-(«(i)>| = q( ^(J)^_^) )- (5.42) 

In the last equality, we used the fact that ff^ are at most of order O(x^) on the real axis, as 
follows from their behavior around ±q. □ 



6 Leading asymptotic behavior of log det [/ + V\ 

In this Section, we prove the result of Theorem 2.1; that is to say, we compute the leading 
asymptotic behavior det[/ + V]^*'^ of det[/ + V] up to o(l) corrections in the x — > +co limit. More 
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precisely, we show that 



logdet[/ + y](°^ = 2 r cUy(i)log'[e-U)] + ^log 

'J n -1- 



-1 



+ 



-1 



6.x dp. 



(2qpUy^ 



A- p 



(6.1) 



This result will be obtained by two different methods based on the integration of equations 
(3.13). The first one, which uses the derivative of the Fredhohn determinant over y, is based 
on the uniformness of the asymptotic expansion for the resolvent for y small enough. It is 
worth mentioning that this way is technically quite involved. The second method deals with the 
derivative of the Fredholm determinant over q. Although we have not been able to provide a full 
rigorous proof for it, we would Uke to draw the reader's attention to this method as it is much 
more direct and simple. 



6.1 The leading asymptotics from the y-derivative method 

Due to Proposition 5.3, the proof of the leading asymptotics of the Fredholm determinant from 
the first equation (3.13), 



R(.A,A), 



(6.2) 



only necessitates the use of R'^^^{A,A) defined in (5.37). Recall that R^^''{A,A) has different 
leading asymptotics in the bulk ] ; ^ [ and near the boundary. Let 5 > be sufficiently small. 
Then 

' Rf{A,A) Ae[q-5;q], 



R^''\A,A) 



where 



R'^^(A,A) Ae[-q + 6;q-6], 
[R^^l(A,A) Ae[-q;-q + 6], 

F{A) 



2in{\+yF{A)) 



\2^;^\oge^{A)-2^x\og 



Kp{A) 



P+ - P(^) 
p(A) - p- 



(6.3) 



r'-^I (A, A) = -v(p {v;x[p- p-]) {iv' log x - 2dx [vlog {p^ - p) - 2dx log Kp] 

+ V{\l/{\ + v) + ^ (1 - v)] + g'} + ixvp'r (y; x[p - p-]) + w'p (v; x[p - p-]) , 



Rf^ (A, A) = -v(p (y; x [p+ - p]) [2v' log x - 2dx [vlog {p - p^) - Idx log Kp\ 

+ ^^y [(A (1 + y) + (A (1 - v)] + g'} + ixvp'T (y; x [p+ - p]) - w'p (y; x [p+ - p]) , 
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HereV^(z) = |logr(z)and we have introduced the shorthand notations 

(f (y; t) = ^{-v,\\ -it) (D (v, 1 ; it) , 

p (y; t) = (5iO) (y, I'it) O (-v, 1; -it) + (5iO) (-v, 1; -it) O (v, 1; /f) , 
T (y; = 1 ; -it) <D (y, 1 ; /f) + (^^O) (-y, 1 ; -/f) <D (y, 1 ; /f) 

+ 0(-y,l;-?Y)(5^(D)(y,l;?Y). 

Moreover, in order to lighten the above expressions and similar ones in the following, we omit 
the expUcit dependence on the argument X of the different functions involved (like y, p, and 
their derivatives V, p', etc.). 

We can now split the integration contour into three parts 

q -q+S q-S q 

J /?^«) (A, ^)j = f jR'^l U, ^) + J y U, ^) + / jR? U, ^) . (6.4) 

-q -q -q+6 q-6 

The bulk integral is carried out straightforwardly. The integrals over the vicinities of the end- 
points are more involved. Consider, for instance, the integration over {-q;-q + 6\ 

Using the asymptotic series for Humbert CHF O (A.9) and the equations (A. 10), (A. 11) we 
get that 

^(a;?)-e'™, 

Aia 



Pia;t) + — {l\ogt - - a) - + a) 

r(l - d;"^" 

T{a;t) + e'"'\\ - 



Y{\-a)Y{l+a) ^ ° ^ ' ^ ' \+t 
2ia 



are uniformly Riemann integrable on R"*" in the sense of the definition of Lemma B.l (See 
Appendix B). Using the integration Lemma B.l as well as the estimates for the integrals of r 
and (p (A.IO), (A.ll), we find 

-q+5 -q+5 

/r d/l e'^ 
/?(«)(i,i)di = - — — -{2y'logx-25^[ylog(;7+-;7)] 
J r(y)r(l-y) 

-q -q 

+v'{il,{v) + il,{-v)) + g'} 

-q+S 

F7TF71 -A2log[x{p-p-)]-ijj{-v)-ilj(y) + —- ; 

J r(y)r(l-y)l^ \+x{p-p-, 

-q+6 

J ^ \xip-p-) + l J 



-q+6 

+ 

-1 



-q+6 

+ ix 

-1 



+ rn ' ^ - ^/r (y_) - ijf (-y_)} + o (1) . (6.5) 

r(l - y_)r(y_) 
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Here the o(l) is with respect to the successive limits x6 — > +00 and 5 — > 0. The two terms 
proportional to log x compensate each other. The remaining part of the first three lines of (6.5) 
is an O (<5) and can thus be dropped. The integral in the last two lines of (6.5) is evaluated thanks 
to the second integration Lemma B.2 (see Appendix B). We get, 



—q+S —q+6 



-dA 



+ 



r(v)r(i-y) 

-1 

{-21og[x(p((5-^)-;7_)] + 2-V^(v_)-V^(-v-)} + o(l) . (6.6) 



r(y_)r(l-y_) 

The integration over [q - d;q] can be treated similarly. The result reads 

q q 



y) 



q-6 



+ 



q-6 



y+e' 



r(y+)r(l-y+) 



{-21og[x(;7+-;7(^-5))] + [2-^/r(y+)-^/r(-y^)]} + o(l) . (6.7) 



So that, 



q q 

J (A, A)dA^J 2in(i%F) ['^^' ^ ' '^J '^^ 



q-S 

. fvlog^ 

J \p-p 

S-q 



+ 



jF \ dA 



1 + jF J in 



yvF IP+-P 
■log' 



q—6 



6-q 



+ 



yv-F- 



2ot(1 +rF_) 

yv+F+ 



+ 



2ot(1 +7F+) 



where we used ^ , , ^ = --.,, 

r(y)r(l-y) 2in{\+yF) 

G-function (2.7), it is not a problem to see that 



in{\ + yF) \p - p. 
[2 log [x {p (6-q)- p-)] -2 + i/i (y_) + i/, (-y_)} 

{2 log [x (p iq-6)- p^)] -2 + 1I, (y+) + if, (-y+)} + o (1) , (6.8) 
yF 



Using the integral representation for the Barnes 



1 1 
J (A, A) ^=dy\J vd^loge- dA + log 



G(l,y+)G(l,y-) 
[xip,-p.)f-''-. 



+ 



q q 

2 J dyvdAlogKpdA-2 J [dAdyv\vlog^y^^ dA . (6.9) 
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Now we should recast the last line as a derivative with respect to 7. We have 



2 J ^dyvdji log Kp - [d^dyv] v log i^JZT^^ 



2v2 / 

k'^ (ci; q) / p+ - Up+- p 



^^^^^ yK''-{-q;q)\ 2qp'+ j \ 2qp'_ 
+ ^ [o-dyVcr log K {crq; q) - crva-dy log k (crq; ^)] 

cr=± 

- 2 I dyd^v 



It remains to apply the identity 



Ou)-i/Ou)v(A) 
2U-JU) 



d/ldju 



1/ 



log /c + y log 



log K + V log 



^ - /I 
^ + /I 



q - A 
q + A 



dA . (6.10) 



d^ 



+ ^ [o-SyVo- log K {a-q; q) - avady log k {aq; q)\ . (6.11) 



Indeed, we have for the r.h.s. of (6. 11) 
RHS 



v+dyV - V dyV+ v-dyV - VdyV- 



q — A 



q-\- A 



+ '■ H — I dA 

1 1 



+ J v{n)dydxv{A)\^ 



A- p + A- p - iO 
1 



+ 



+ J yOu)Vw{ 



jO q - crA - iO 



1 1 

+ 



dAdp 



dA 



iA-p + iOf (A-p- iOy 



dAdp. (6.12) 



There we have regularized all the integrands and then performed an integration by parts. On the 
other hand, one has for the l.h.s. of (6.11) 



1 

LHS = J [dyv (p) dxv {A) + v(p) dyd^v (A)] | 



1 



iO A-p - iO 



dAdp. (6.13) 



Taking the last integral by parts we arrive at (6.12). 
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Thus, the l.h.s. of (6.8) is presented as a derivative with respect to y. Since the asymptotic 
expansion is uniform in y we can integrate this result from to y. As log det [/ + V\ |y=o= we 
get the desired result. 



6.2 The leading asymptodcs from the ^-derivative method 

The method we use here is based on the second equation in (3.13), 

dq log det [/ + y] - /? {q, q)+R {-q, -q) . (6. 14) 

For the purpose of this sub-section, we assume that |2l (v(/l))| < 1/4. Indeed we are then able to 
use the pointwise estimates for the resolvent estabUshed in Proposition 5.5. Such a restriction on 
|!R (y {X))\ could be relaxed by much more refined estimates. Recall that one has for A uniformly 
away from the boundary Sn corresponding to disks of radius 5, 



\R {A, X) - R^^^ (A, A)\ < < , with e = 2sup l^lvl , 



(6.15) 



so that the q anti-derivative of R (q, q) +R (-q, -q) - R^^^ (q, q) - R^^^ (-q, -q) will be a o (1) in 
the X +00 limit. 

Equation (6.3) allows us to determine the value of 

7?(0)(^,^):._(F«0)(^)|(^,F«^(«)(i)) 



= ^fiyA)f!:'\A) {5,logA -5,log/_} 

at both endpoints q and -q. 

Consider, for instance, R^^^ (-q, -q). We have, for A e D-q^s, 



(6.16) 



{A, A) = -v{A) O (-y, 1 ; -i^.q) O (y, 1 ; i^.q) |25^ log e^{A) 

- 25a[ log Kp{A) + v{A) log ^q] + v'(i)[«A (1 + y) + (A (1 - v) ] 
" (5,(D) (y, 1 ; /^_,) (5,(D) (-y, 1 ; -i^-q) ' 



- ixp'(A) 



-v'iA) 



(D(y,l;/^_,) cD(-y,l;-/^.,) 
(5i(D) (y, 1; /^_,) (5iO) (-y, 1; -/^-,) 
a)(y,l;j^_J a)(-y,l;-j^_J 



(6.17) 



where, so as to lighten the formula, we have omitted the argument A of v(A) when y appears 
as an argument of another function (here (A or O). The symbol 5, stands for the derivative of 
a CHF with respect to its variable, whereas di stands for the derivative with respect to its first 
argument. Recall also that ^-q - x[p(A) - p (-q)] and (q = x[p (q) - p (A)]. 

It is remarkable that the last two terms involving derivatives of CHF vanish in the A ^ —q 
limit. The resulting expression can be further simpUfied thanks to the identities: 



log Kp(A)^ log k(A) + v(A) 



{.og(. 



piA) 



log 



A + q 



p{A) - p 



:)}• 



v{A) v'{A) [il,{l + y) + - y)] = ^;^ log G(l, y) + 2v{A) v\A), 



(6.18) 
(6.19) 
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Thus, we obtain 

R^'>\-q,-q) = -2y_[5^1oge+U)]| + 2yly_logx 

+ 2v'_v- log(2qp'_) - — + vl^- 2y_y'_ + 2y_ [5^ log k{A)\ i , (6.20) 

q P- \x=-q 

where we have used the notations (2.3), (2.4). 

The final aim is to integrate (6. 14) over the variable q. One should keep in mind that the 
function k{A) = k{A; q) is actually a function of the two parameters X and q. Therefore, one 
should replace partial X derivatives dX.A-±q by total q derivatives thanks to 



— [log K{-q;q)\^ -dx log k{A; q)\ + dq log k{A; q) i 

U=-? I 



Then R^'^^ (-q, -q) is almost a total q derivative: 



R^^^ i-q, -q) = -2y_ [dx log e^{A)] , + log 



G(l,y_) 



(6.21) 



(2qp'_xr- 



2y A [\ogK{-q;q)\ + v ^^ "^"^ (6.22) 
aq q 



Similar calculations based on the expressions (5.20) for ff^ around q lead to 



{q, q) = -2y+ [d^ log e+{A)] , + log 



<:?(!, V+) 

(2^/^;x)^' 



+ 2y+ [log ^ (^; q)] - y+ ^^^^-Jll . (6.23) 
d^ q 



Hence, we have 

dg log det [/ + y] = 2 y y^[dx log ^_(^)] | + T" 



G(l,y+)G(l,y„) 



(TVn 



{2qp'_xY- (Iqp'+x) 
d\og K(crq;q) (y+ - y.)^ 



dq 



+ o(l). (6.24) 



It remains to express the last Une as a total ^-derivative thanks to Lemma B.3. After an 
integration with respect to q of (6.24) we arrive to 



log det [/ + y] = 2 J dAv(A) log'[e_(^)] + log 



G(l,y.,)G(l,y-)^'--(g;g) 
{2qp'^xy'i2qp'_xf-KH-q;q) 

1 f v'{A)v{p)-v{A)v'{p) , ^ , 

+ - I dAdu +C + o(l) , (6.25) 

2 J A- p 
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where C is a ^-independent integration constant still to be determined. 

One can give arguments that this constant should be also y-independent. Indeed, the asymp- 
totic expansion of the Fredholm determinant, being a functional of the holomorphic function 
yFiA) in U , can depend on this function either in the integral form with integration over [-q; q], 
or through the values of -yF and of its derivatives at the ends of the integration contour -q and 
q. In both cases the result should depend on q. Hence, the (jr-independent constant C can not 
depend on yFiA) and, thus, it is y-independent. We can then fix the constant C by setting y - 
in the asymptotic formula. This yields C -0. 



6.3 The first corrections to tlie leading asymptotics of tlie Fredliolm determinant 

We close this section by deriving the sub-leading corrections from the x-derivative (3.19) of 
log det [/ + V]. This will constitute the proof of Proposition 2.1. 

In order to prove the claim of the Proposition 2.1, one has to derive the first two sub-leading 
corrections for the matrix IT. As one might expect the computations are, by far, simpler that 
those necessary to fix the constant. We also would Uke to point out that one can obtain the 
sub-leading asymptotic of det [/ + V] by the (^-derivative method. However, the computations 
are quite involved, so we omit the presentation of this method. 

We derive the first term in the l/x expansion of log det [/ -i- V] thanks to (3. 12) : 

log det [I + V]^ ^ ^ p{A) tr [5^ (A) o-^x'' (A)] , (6.26) 
r([ -?;?]) 

where we chose the contour T ([ ; ^ ]) to lie outside of the contour En but still in the region 
of holomorphy for p. There the solution for the RHP for has a simple form: 

X(A)^U(A)a-''HA). (6.27) 

In order to derive the first correction to the leading asymptotics, it is enough to consider the first 
two terms in the asymptotic expansion for n(/l): 

nW) = /,.S!!l»,Sf^,0(.»-»). (6.28) 

JC X 

There, as follows from (5. 15), the O is uniform on the whole contour T ([-q;q]). Thus 

log det [/ + V] = - ^ g p (A) ^ + 1 ^ ^ p (^) tr [cr3 ^.H^^) (A)] 
r([ -?;<?]) r([-^;^]) 

+ ^ ^ ^ piA) tr {^3 [d,n^'\A) - n«(i)5,n«(i)]} + O (x^^^-^)) . (6.29) 
r([-9;9]) 

The first term in this expansion will yield the leading correction. Indeed, 

g 
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Here we shrunk the contour to [ ; ^ ] and used the jump condition for or. 

In order to evaluate the higher order corrections in (6.29) we need to derive the expressions 
for the matrices U^^^ = Csn[A^^^] and n^^^^ = C^nLn+^A^i) + A^^)] outside of Sn. An elementary 
computation of residues yields: 

5^n(i>(^) = — — , (6.31) 

as well as 

5^n(2^ {A) - n^i) {A) (9^n(i) {X) 

(dA^fl^ {(rq; X) + A;^^^ (crq; x) (s^aJ^J {aq; x) 



cr 



[(A - a-q) p'a-f 
2^f^^{cTq■x) + [^^^^{cTq■x)f ( p'; 



2{p'J \pU^-(rqf {A-aqf 



[^^^^{q■x),^^^^{-q■x)\ 



(6.32) 



2qp'^p'_{A-q){A + q) 
Thus the 1 jx term in (6.29) gives the coefficient in front of log x in (2.6). Indeed, 

^ ^ p (i) tr [era ^^n^^^ {A)] = ^tr [craAjl^^ {-q; x) + (t^a'^I^^ {q; x)] 

^-(vi + vl). (6.33) 

We now focus on the last term in (6.29). It yields, after an x integration, the first correction 
to (2.6). A straightforward computation leads to 



|i p (A) tr (era [5,n(2) (A) - U^'^ (A) ^;^U^'^ (A)]} 



§iin 
r([-?;?]) 

P+-P 



o-=± 

The first term corresponds to the oscillating correction: 

tr {era [a[|5 (,; x) , A^^^ (-,; x)]) = 2v.v_ - ^i^j . (6.35) 



The last term gives the non-oscillating one: 

dy^ 

tr [a^ {d^f^^) {crq; x)} = -2-^, (6.36) 
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and 

Z 77^^ h ^) (^-^^S) 

o-=± P"- 

2 2 
= — v^5a (log u {A; x))\ + —vld^, (log u (A; x))\ 

P+ U=9 P- \A=-q 

2^{w,log. + 4log,<. + ^^A(4)_r^| , (6.37) 

hi P-r [ V^d^\P^r/ P- J 



cr=± 

with 

= {^f, (6.38) 

«- = ^'^"'^ {(2^y_r^(-^; q)}-^ . (6.39) 

Putting all this together we obtain 

^.logdet [/ . V] = -/ f di v(i)/ iA)-i^- 

^^^/'^r I d^ d^\K/ ? J 

(2qfp'p'_x^ \u- J 



ijx(p+-p-) 1 

1-^,^1. (6.40) 



The first two terms reproduce the already known answer for the leading asymptotics. The re- 
maining ones reproduce the first oscillating and non-oscillating corrections as given in Proposi- 
tion 2.1. Note that for the oscillating terms, one only should integrate the exponent with respect 
to X as all the other terms will give subdominant contributions. 

7 Applications to truncated Wiener-Hopf operators 

Truncated Wiener-Hopf operators appear in many domains of mathematical physics such as 
scattering or diffusion processes. Moreover, many observables (dressed energy, momentum 
or dressed charge) related to quantum integrable models are solutions of integral equations of 
truncated Wiener-Hopf type (2.17). 

Let us recall that a truncated Wiener-Hopf operator can be interpreted as an integral operator 
/ -I- A" on L-^ (R) such that it acts on !?■ (R) functions according to 

X 

(I + K).(p(t)^(p(t) + J dt'K{t-t')<p(t')dt'. (7.1) 
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The kernel K is characterized in terms of its Fourier transform F: 

m = r-' it) , with r-'{Fm = ^ J d^^(^) e-''^ e L'm. g.i) 

R 

The study of truncated Wiener-Hopf operators is equivalent to a 2 X 2 matrix RHP. Another 
facet of this equivalence is the correspondence between a truncated Wiener-Hopf operator and 
the GSK acting on R. in which p - id and g - 0. Indeed, it is easy to see that 

K.<p^r-^ oVor[(p], (7.3) 

where V acts in (R) with a kernel 
The operator identity 

I + K^r-\(I+V)r, (7.5) 

together with the facts that V is trace-class and are continuous, ensures the equality between 
the Fredholm determinants: 

det [I + K]^ det [l + v]- (7.6) 
The kernel V is related to 

T/^t ^ l irr^Mrr , siu (x (g - ?/) /2) 

V (^, rj) = y^F (^) F {?]) — (7.7) 

by a similarity transformation. Hence, 

det [/ + K]l2(o,,) = det [/ + V]^2(r) . (7.8) 

7.1 The Akhiezer-Kac formula 

Our study of the generaUzed sine kernel allows us to recover the Akhiezer-Kac formula de- 
scribing the large x behavior of Fredholm determinants of truncated Wiener-Hopf operators: 

Theorem 7.1 (Akhiezer-Kac [1, 30]). Let I + K be a truncated Wiener-Hopf operator as above 
and such that 

• F is analytic in an open neighborhood U ofR; 

• F goes sufficiently fast to at ±oo ; 

• \ + F {^) does not vanish on U. 
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Then 

oo 

log det [I + K]^ xt(0) + E[F] + o(l), with E[F] = ^ tT{t) T{-f) At, (7.9) 



in which 

T(t) log(F(^) + 1) e-'^^ . (7. 10) 

R 

Proof — The large x asymptotics of det [/ + K] follows from (7.8) after taking the q +oo 
limit in the leading asymptotics for the corresponding generalized sine kernel (2.6). This limit 
may seem a little heuristic as we did not specify any estimates in q for the small o terms with 
respect to the leading asymptotics. However, the validity of such a limit may either be seen 
by refining all the estimates obtained in the previous section or by considering the RHP for x 
(3.1) on the whole real line from the very beginning. We shall make the second approach more 
explicit in the forthcoming subsection 7.2. Here we formally take the q +oo limit in the 
leading asymptotics of Theorem 2.1. 

One should notice that, in the asymptotic formula (2.6), all the terms evaluated at the end- 
points vanish due to the fact that y(±^7)log^ — > when q — > +oo, which is a consequence of 
the sufficiently fast decrease of F at infinity. Hence, the only constant contribution E[F] to the 
asymptotics of log det [/ + K] is given by the integral 

Ein^ lim i f-^WW-'W'Wd.d^. (7.11) 

q->+oo 2 J A- jX 

-1 

Let us recast the constant E [F] in a more standard form. We have 

log' (F(^) + 1) log (F(ij) + 1) - log' (F(77) + 1) log (F(^) + 1) 



E[F] = y d^dT] 



d^ drj dx dy { i + i 1 T{x)T{y) {x-y) e'^+'^^ 



Let H be the Heaviside function, then 
E 



[F] = - ^ dj]dxdyT (x) t {y) {x - y) (e'''^^+^^// iy) - e'''(^+y>// (-y)) 

R 

+00 

= -i/d,rWT(-,)(-2,).l/d,rWr(-,)(-2,) 





+00 



= ^ dyT{y)T{-y)y , 
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which ends the proof of Theorem 7.1. □ 

It happens that this correspondence between truncated Wiener-Hopf operators and gener- 
alized sine kernels can be pushed further so as to obtain the asymptotic behavior of Fredholm 
determinants of truncated Wiener-Hopf operators with symbols having Fischer-Hartwig type 
discontinuities. Considering the GSK for finite q somehow corresponds to the asymptotic be- 
havior of a determinant whose symbol has two jumps. The case of symbols having general 
Fischer-Hartwig type singularities is studied in [33]. 

7.2 The resolvent of truncated Wiener-Hopf operators 
Proposition 7.1. Let I + K be a truncated Wiener-Hopf operator on ] -x ; x [ , 

X 

[{I + K).g\{t)^g{t) + J K(t-t')g{t')dt', with K(t)^r-'[Fm. (7.12) 

-X 

Suppose that there exists 6 > such that 

• F admits an analytic continuation to |z : |3(z)| < <5|; 

• ^ F ± i6) e (W); 

• the analytic continuation ofl+F does not vanish on U. 
Then the resolvent I - R of I + K fulfills 

R^,)^ f ^F(^)|^e'-^(^-)-^e-(^-)|^.0(e-n (7.13) 



R 

where 

ail) = exp J r ^^d/i 1 , and v{A) = — log (1 + F{A)) . (7.14) 
\J li-A \ 2n 



Proof — The GSK associated to I + K through the transformation T'^ olI + Vlof^-I + K 
acts on the whole real axis with the kernel 

J(^-r])x _ J(r]-Ox 

V (^, 77) - F{^) —r— — . (7.15) 

Just as for the leading asymptotics of log det [/ + K\ (see Section 7.2), one can obtain the 
leading asymptotic of the resolvent of W just by taking formally the limit q — > -i-cx3 in all the 
expressions derived in the former part of the article. Note that in this process all power law 
corrections vanish: they are computed as contour integrals around ±q and, since F approaches 
sufficiently fast at infinity, the residues at ±q vanish in the q +00 limit. However, in order 
to justify this hmit, one should also check that all the uniform estimates still hold for q — > +00. 
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An alternative way is to consider from the very beginning a RHP for x on the whole real 
axis R. This is actually much simpler, than the RHP on a finite interval. Then it is enough 
to perform the first two transformations described in Section 4 so as to obtain jump matrices 
that are already uniformly close to h up to exponentially small corrections in x. Moreover, the 
jump matrices for this RHP are given by M+ and MZ^ (4.11), so that they approach the identity 
matrix at A ^ oo just as fast as F goes to zero at infinity. As expected, there is no need for 
parametrices, and the corrections are immediately exponentially decreasing with x. It means 
that, up to uniformly exponentially small corrections, the resolvent Ry of V is given by 

< ^) - i"-^/'''"' - ^e-^C^-)} , (7.16) 

where as usual a {A) is given by (7.14). Note that the integral in (7.14) is well defined in virtue 
of the assumptions made on F. 

We should now take the Fourier/inverse Fourier of Ry in order to get R. To this end, we must 
justify that the sub-leading corrections do admit a Fourier transform in two variables. Recall the 
exact expression for the resolvent: 

Ry iA,f^) = iA,f^) + < (A) lElS^lIMzIliF'^m (^) ) . (7.17) 
^ A- n 

Thereby, the matrix n is defined in terms of n+ {A), the limiting value of H on Sn when A 
approaches a point of £n from the "+" side of the contour: 

'(Z) 



+ 



The integrabihty of F as well as the asymptotic condition H {A) — > h guarantee that the 

/I— »+oo 

integrals are well defined. Thus one readily infers from (7.18) the asymptotics of H on the real 
axis: 



n(i) = /2 + — ^ + 0^— j , (7.19) 

where C is some constant matrix and where we have explicitly extracted the exponential decay 
in X of the matrix C. Hence using the boundedness of f^^ on the real axis we obtain that 



Ry (A, IX) = {A, IX) + ^ '"C^ + o ^ j-^j . (7.20) 

Hence the corrections admit a Fourier transform in A and an inverse Fourier transform in ju as 
oscillatory integrals. Therefore, taking the Fourier transform does not change the nature of the 
corrections. □ 
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It is clear that, up to a similarity transformation, a Wiener-Hopf operator on[-a ;b] has the 
same generalized sine kernel than the same operator acting on [ ; a + ^ ]. Therefore our method 
works for any interval, of course up to a similarity transformation on the resolvent (7.16) of V. 
We chose here to present this less standard form of Wiener-Hopf operators as it fits better the 
forthcoming application. 

We apply our asymptotic inversion formula for truncated Wiener-Hopf operators acting 
on a symmetric interval [ -;c ; x ] to re-derive some formulas concerning the low magnetic field 
behavior of the so-called dressed charge [6]. This function, traditionally denoted Z{A), describes 
the intrinsic magnetic moment of an elementary excitation above the ground state in the XXZ 
spin- 1/2 model. It satisfies the integral equation: 

X 

Z{X)+ {ahK{X-h)ZUi)^\, with TO- ^ -w^r-rww, (7-21) 
J 2;r smh(/l + ?^) smh(/l - iQ 

-X 

Thereby K is often called the Lieb kernel and ^ e ] ; tt [ is some real parameter describing the 
coupling of the model. The large parameter .jc is a function of the external longitudinal magnetic 
field; it goes to infinity when the magnetic field vanishes. 

For the study of Z, one should distinguish two domains in the interval [ -x ; x]: the bulk, i.e. 
the region |/l| « x, and the boundaries X ~ ±x. While the asymptotic value of Z{X) in the bulk 
(|/l| «; x) is enough to describe the intrinsic magnetic moment of elementary excitations, the 
value of Z at the boundaries (±x) determines the critical exponents of the two-point functions 
of the model [20, 21, 22, 5]. As we will see, the bulk and the boundary behavior of the dressed 
charge differ fundamentally. 

First, let us note that, setting directly x - +oo in (7.21), one can solve explicitly the integral 
equation for Z by taking the Fourier transform: one obtains in this case that Z{A) is equal to a 
constant value Z{X) - nj [2 (n - ^)] on the whole real axis. 

Let us now consider the limit x ^ oo in (7.21) in a more accurate way, namely, taking x 
large but finite, and use the method described above. Let K be the Fourier transform of K, 

ftf)^rmg)= ^'°".'l<f-f" . (7.22) 
smh(^^/2) 

Then in virtue of Proposition 7.1, 

X 

Z(A)^l- J R(A, n) dn (7.23) 

-X 

J 2m ^ [a-(0 a_(0) J 

R 

First let us study the bulk hmit i.e. \A\ «: x. Using the jump equation satisfied by a±: 
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[1 + K(A)]a+(A) - a-(A), we recast the integrand as 

Z(^) . 1 - r d£^3f)^_ _ 

J 2/7r[i +^:(^)](^_ jO+) a^-(O) 

i(/2) a+{0) ei^-m-^m 



1^ r Kj^+ig: 



•(0) J 2ot 1 + + a+{^ + /^/2) ^ + 
1 + ^(^ - j^/2) «-(^ - '^/2) ^ - il2 



(7.25) 



Here we have separated the integral into two parts and then moved the contour to the upper/lower 
half -plane. This gives a pole contribution from^ = /O. The integral appearing in (7.25) is clearly 
a O [Q-i^-mi^"^^ So that, in the bulk, 

Z{X) = — ^ , (7.26) 

1+^(0) 2(7r-^) 

up to exponentially small corrections. As expected, we recover the value of Z obtained in the 
case of an infinite interval. Note that the corrections become larger and larger as we approach 

any of the endpoints ±x. 

Let us now study the behavior of the dressed charge at the boundaries. Since the kernel K 
is even, so is Z. We can thus focus on a single boundary, say A = x. We have, 

Z(.) = I - f \m J^'^ M . (7.27) 

J 2j7r(^-?0+)l a_(^) \+K{^)a-{0) I 

R 

As before, the integral of the second term gives an exponentially small contribution O 

The integral of the first term is expUcitly computable. Using once again the jump equation 

satisfied by a±{^), we have, 

J 2in ^ - iO ^ > 

R 

= a+(0) + O (e-^^) . (7.28) 

We have computed the remaining integral by residues, since a'^{£;)-\ = O for ^ ^ cx5 in 
the respective half plane of holomorphy. 

For an even kernel like the Lieb one K{^) - K{-^), and then it follows from the integral 
representation (7.14) of a that q;+(^) = aZ^{-^). This means that 1 + K{Qi) - a+'^(0). Hence, for 
X large enough 



Z{x) ~ ^ - -, (7.29) 

Vl+/s:(0) V2(7r-^) 
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and the value of Z(A) at the boundary is the square root of its value in the bulk up to exponentially 
small corrections. In the Umit x — > +00 this correspondence becomes exact. 



8 Asymptotics of multiple integrals 

We have already mentioned that the asymptotic expansion of the Fredholm determinant of the 
GSK can be used for the asymptotic analysis of correlation functions of quantum integrable 
models. For a relatively wide class of integrable systems, the correlation functions can be 
presented as series of multiple integrals of a special type [31]. These series can be summed up 
to Fredholm determinants for the models equivalent to free fermions. In the general case, such a 
reduction to determinants is not known. However, the asymptotic behavior of individual terms 
of the series can be derived from the asymptotics of the Fredhohn determinant of the GSK. In 
the present section we consider this problem. 

More precisely, our purpose is to derive the large x asymptotic behavior of the following 
type of integrals (cycle integrals): 

Inirn] = (p — Tn\ \ ' 7 ^, r • (8.1) 

In this expression, Tn is a holomorphic function of 2n variables /li , . . . , /l„ , zi , . . . , z„ in J7" x W, 
in which U and W are open neighborhoods of [-^;^], and r([-^;^]) denotes some closed 
contour in W surrounding [ ; ^ ] with index 1 . We moreover assume that Tn is symmetric 
separately in the n variables Ai,...,Xn and in the n variables zi, . . . , Zn- Finally, we agree upon 

8.1 Leading asymptotic behavior of J„ \Tn\ 

Let us first suppose that the function Tn is of the special (factorized) type 

j^j )-n[^(A,)^(z,)], (8.2) 

^ ' (=1 

where ^ is a one-variable holomorphic function in U and is a one-variable holomorphic func- 
tion in W, non-vanishing on W. For two such functions (p and (p we introduce the associated 
GSK given by (1.6) provided the identification F{X) = (p{X)<p{X) and e^^^^ = ^(z) is made. 
Then, the integral (8.1) can be expressed in terms log det [/ -i- as 



= ^ log det [/ + V^f'^^] I . (8.3) 

(n-1)! ^ |y=o 

In this specific case, it is straightforward to obtain the leading asymptotic behavior of the mul- 
tiple integral (8.3) in the large x limit thanks to the results of the previous sections. 
This remark leads us to the following definition: 
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Definition 8.1. Let U, W be two open neighborhoods of[-q;q \ and let 9i{U) ( resp. ^{W) ) 
be the set of holomorphic functions on U (resp. on W). Let also 

^n'^ = I X ^i""'^'^ -^P^^, (.<P{, H) e X "KCW) and(l>e\ ^ ol , (8.4) 

in which ^^^'^^ denotes a pure factor function of2n variables defined in terms of (p and (p as in 
(8.2). We define the linear functional I^^ on S^'^ as 

/f [T^^'^'] = a; log det [/ + V^-^'^Y") I , (8.5) 

(n-1)! |y=o 

and by imposing linearity on functions <j^j^'^''^f\ Here V^'^'^^ denotes the generalized sine 
kernel (1.6) with F{A) = (p{A)(f>(A) and e^(^> - (/>(z), and log det [/ + y^*^''^)]^^^ denotes the leading 
asymptotics of the Fredholm determinant log det [/ + V^'^'^^j as in Theorem 2.L 

It is easy, using the expression (2.6) of log det [/ + V^'f''l'^f\ to obtain an explicit expression 

= ^ ^ip\A)r-\x){ixp'{A)m+4>'m 

^ {bn - Cn log {2qp'^x) ) Vp{o-q) <p{o-q)f 

[ip((Tq)4>{crq)r - [<p{crq) 4>{crq)]P [ip{A) cpiAJf-P 



+ 

cr=± 

n-1 

' dA 



^ II— ± n 



^=+„=i- pin- p)(q-o-A) 

I -a 



n-l ^ 



where b„ and c„ are given by (2.23). The x — > +00 asymptotics of the Fredholm determinant are 
uniform in 7 to any fixed order n in 5". This means that 

In[rn-I^n\rn+0{\). (8.7) 

In the next proposition we show that 7^*'^ can be extended into a Unear functional on the space of 
holomorphic functions T„ (not necessarily of the form Tn^''^^) that are symmetric in n variables 
A\,. . . ,An and n variables zi, - ■ ■ ,Zn separately. This extension of 7^*^^ , as we prove below, is the 
good way to evaluate cycle integrals (8.1) with such arbitrary symmetric functions 'F'n. 

Proposition 8.1. Let U and W be open neighborhoods ofy—q\q], and let SyvciJJJ, W) be the 
set of holomorphic functions Tn on [/" x W" of In variables Ai, . . . ,A„,zi, ■ ■ ■ ,Zn, symmetric 
in the n variables Ai,...,An and in the n variables z\,...,Zn separately. Then I^^ extends 
to a continuous linear functional on Syrn^iU, W) endowed with the topology of the sup norm 
convergence on compact sets. 
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Proof — contains at most first order derivatives of the functions (p and 4>- Now recall 

o 

that, for any compacts K, P such that^ K <z P and P <zU 

\fk e N, 3cfc such that V0 e <?/(f/) , \\<p^%.j^ < Ck ||0|Io;p , (8.8) 

where \\.\\[)-k = sup^.^^ |.| is the sup norm with support on the compact K. In consequence, 7^*''' is 
continuous on S^'^ . The latter is dense in S^'^, with 

S^'"^ = |g Ti""'^'^ ■,peN, (cpe, cf,e) e 9<iU) x <H(Wo| . (8.9) 

Hence 7,,*^^ extends by density to a continuous linear functional on <S,f'^. Due to the density 
Theorem C.l (See Appendix C), we have that S^'^ is dense in ,Sym„(t/, W). Therefore 7,^°^ 
extends to a linear functional on Sym^{U, W). □ 

Corollary 8.1. Let U and W be open neighborhoods ofY-q;q\, and let Tn ^ Syra^{U, W). 
Then, 



+ {bn - Cn log {Iqp'^x}) Tn ( j^^jl 



e=0 



n-l % Jn\ f in ] - 



i'^^fhild'J p{n-p){q-(rX) 



p 

n-l ^ 



(^ 



2 (2;r)2 ^ j (n - p) (^ - //) j^^^" \{A + e], {A]P-' , {/z}""^ 



^.1 (n- 



i/. + 6i,i;.ri,{^r/' 



e=0 



There {/I}" denotes the set formed by n copies of the same parameter X. 



Proof — Apply Theorem C.l to (8.6). □ 

Finally, we have the following large x asymptotic behavior for integrals of the form (8.1) 
(which seems hardly attainable through a direct analysis of the multiple integrals): 

Theorem 8.1. Let U and W be open neighborhoods of \ —q\q\ and let 'Fn ^ Syva^{U,W). 

Then, when x — > +cx5, the integral In{%i\ (8.1) behaves as 

J„r„]=7f[r„] + o(i^), (8.11) 



Here P is the interior of P 
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the explicit expression ofI^\'F„] being given in Corollary 8.1. 

The whole difficulty of the proof is to show that the small o (1) in (2.5) is preserved by the 
density procedure formulated in Theorem C.l. This is nontrivial since the series converging to 
'J^n may not converge absolutely. We need therefore, so as to prove this theorem, to study more 
precisely the sub-leading corrections and to see how they pass through all the steps described 
above. This will be done in the next subsection. 

8.2 Study of sub-leading corrections 

In this subsection, we study the behavior of the sub-leading corrections to logdet[/ + V]^^^ when 
the above procedure is applied. In particular, we will show that they indeed remain subleading, 
which will prove Theorem 8.1. In fact we will prove an even more general result: 

Theorem 8.2. Let U and W be open neighborhoods of[—q\q \ and let T'n e SymJJJ, W). For 
any positive integer M, there exists a continuous linear functional I'i^^ such that 

Jnr„]=/rr„J + o|^j when;c^+oo. (8.12) 

The explicit expression for 1^ can be obtained by some perturbative computations that 

become more and more involved with the growth of M. We will nevertheless obtain the general 
structure for l''n^\ showing that it can be decomposed in terms of non-oscillating and oscillating 
contributions, with oscillating factors of the form m e Z*: 

N=\ meZ* N=2\m\ ^ 

\m\<M/2 

j(N, nose) i^^j (resp. YF'nT) being given in terms of the function T'n and of its derivatives 

up to order (resp. up to order N -2) evaluated at ±q or integrated from -q to q. 

8.2.1 General strategy 

In the previous subsection, we have defined the functional ij^^ from the leading asymptotic 
part logdet[/ -i- V]^^^ (2.6) of the GSK. More precisely, we have seen in Corollary 8.1 that 
a^logdet [/ + V]^^^ \y=o yields the functional (-1)""^ (n - 1)! ^ [fn] after the density proce- 
dure, as explained in Proposition 8.1 and Theorem C.l, is applied. In order to estimate the 
corrections to lf'\'F„] for the large x behavior of cycle integrals J„[!F„] of length n (8.1), we 
have to take into account the corrections log det [/ -i- V]^^^ to log det[/ -i- V]^^\ 

log det [I + V]= log det [/ + y]*°* + log det [/ + V]'"'' , (8.15) 
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and to analyze the effect of the density procedure on the n-th y-derivative of the subleading part 
dy log det [/ + Vf"^ \y=o. We will show in particular that it preserves the small o(l) with respect 
to the X +00 limit, i.e. that d" log det [/ + Vf"^ |y=o can only generate o(l) corrections. 
In the spirit of Definition 8.1, we therefore introduce the: 

Definition 8.2. Let U, W be two open neighborhoods ofY-q;q\. We define the linear functional 
If" on S^'"^ as 

/r'ri^'^^] = log det [/ + V^'^'^T^ , (8.16) 

{n - 1)! |y=o 

and by imposing linearity on fiinctions H^^j (f^(V(4e)^ Here, as in Definition 8.1, ^^^''1'^ denotes 
afactorized fiinction ofln variables defined in terms of (p and <p as in (8.2), and V^f-I''^ denotes 
the generalized sine kernel (1.6) with F{X) = (p{X)<p{A) and ^^"^ = (piz). 

According to the scheme presented in the previous subsection, the next steps will be: 

• to obtain a convenient representation for P^^\T'n^''^\ this means in particular to obtain 
the form of n-th y-derivatives of log det[/ + V]^"'' in terms of the functions F and g, to set 
g{z) - log (p{z) and F{z) - (p(z) (p(z), and to estimate this result in the large x limit; 

• to apply the density procedure: one should first extend by density and continuity the 
functional 7^"'' to the space S^'^; then, for any holomorphic function fn in 2n variables 
Ai,...,An,zi,...,Zn, symmetric separately in the variables A and in the variables z, one 
has to consider a sequence {(pk, (pu) such that Yjk=\ Tn^'"^''^ Tn so as to be able to define 
and characterize and to see how it behaves in the large x Umit; 

• to refine the procedure in order to get an asymptotic expansion of l^n^\fni- 
S.2.2 7-derivatives of log det[/ + V]'"'' 

As in Section 6.3, we will obtain the corrections to logdet[/ + V]^^^ through the x-derivative 
path, starting from formula (3.12) that we recast as 

1 

log det [/ + V] = -i J dA p'(A) viA)+ ^ ^ p(A) tr [[d^nU)] as U'^A)] . (8.17) 

Here, as in Section 6.3, we have chosen the contour T {[-q ■,q]) to lie in U and to encircle Sn, 
which means that;^f(/l) - Il{A) a~"'^ (A) on T {[-q;q ]). Integrating this equation with respect to 
X, we obtain 

X 

log det [/ + Vr"" = Jdx' p{A) jtr { [^^n (A)] cr^ H"! (A) ] 



+°° r([-q;q]) 

.if!! 

2in(A-z) 



X' J '^:-^-> 
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The convergence of this integral will be proved later on. We recall that the second term in 
(8.17) produces also, when integrated over x, the logx term appearing in the definition (2.6) 
of logdet [/ + V]*^"^ (see (6.33)). We have therefore substracted the corresponding contribution 
(second term of (8.18)) in the definition of log det [/ + V]™''. 

In order to obtain the n-th y derivatives of this expression, hence we have to compute the 
y-derivatives of d;J\{X) and of Yl~^{X), which in their turn follow from those of the jump matrix 

• 7-derivatives of A 

In order to determine the n-th y-derivative at y = of the jump matrix A (z), it is convenient 
to express it in the following form: 



(8.19) 



Here, the matrix A(z) depends on y only through the combination y^"(z), whereas k-""^ (z) de- 
pends on y through the combination da [y (z) - vUjl)] /{z - u). 

It is easy to compute the multiple y-derivative of A(z) at y = 0. It is given as 



d; A(z) I = Ad,.3 ,«/2 [dlAoiz)] I • F\z). 
lr=o lr=o 



(8.20) 



In this expression, AdxLJ'] stands for the usual adjoint action of the matrix X on the matrix Y, 
and Aq denotes the jump matrix A at F = 1 and g = 0. 

It remains to compute the y-derivatives of a:*""' (z). They follow from the Faa-di-Bruno 
formula: 



i>'-o pu-,Pn=o rip.! .=1 



(-1)-^ rF'iz)-F'(p) 



2ins 



f 



z-ju 



(8.21) 



S^=iiPs=n s=l 

Therefore, gathering these informations and applying Leibnitz's rule, we obtain that 
^rAW, = E E Z 

Ir-U p+po+q=n pi,...,Pn=0 qu...,q„=0 (pj! (q,)\ 
P0>1 ^^^sp,=p J.1^isq,=q s=l 



X Ad,.3,fe,/2 [ (-0-3)^-1^^ • 5f Ao(z)| • (crs)^-!"' 
L 7=0 



FP\z) f] 



-\Ps+qs 



{-ly rF^iz)-F'(ju) 



= J] 6(f^'%;x)F^''(z) f] 



po+'^sps=n 
Po>l 



s=l 



lins 

q 



z-n 



-q 



F'{z)-F'{^) 
z-n 



Ps 



(8.22) 
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Note that, in the first line, we have extended for convenience the sum over parameters and qs 
up to n, since anyway, due to the constraint, Ps < p and ps - if s > p (resp. qs < q and qs -0 
if s> q). 

In the last line, we have changed the order of summations and incorporated all the F inde- 
pendent prefactors into the definition of the matrix 8^*''''^ (z; x). More precisely. 



s=l 



{-If 

lins 



S\Ps 



p+q=n-po S^^ii9,=3 n (qsV- iPs " qsV- 
s=l 



xAd 



igo-a «(z)/2 



Is 



lr=o 



. (8.23) 



From the properties of the jump matrix Ao(z), it is easy to see that the diagonal entries of the 
matrices 8^''''') (z; x) are a 0{x~^), whereas their off-diagonal ones are a O (log^" xlx) uniformly 
on the contours Xn- 



• y-derivatives of 

Let us recall the integral representation for 5^11 (/I), which is a direct consequence of (5.5), 



(8.24) 



Therefore, the y-derivatives of d^Yl{X) can be directly obtained from the ones of A(/l) and of 

n+(i). 

Recall that 11+ (/I) satisfies the following integral equation on M2 {l? (£n)): 



where the operator C^^ is defined by 

[M] = C^^ [MA] , VM e M2 (l^ (Sn)) . 



(8.25) 



(8.26) 



This matrix Cauchy operator is invertible, at least for x large enough. Indeed, using the conti- 
nuity of the scalar Cauchy operator: 



3 C2 > such that, Vg e (Sn) , ||C^^ ys\\\j2^^^) < ci Ml^^^) , 

fulfills: 



one gets that the operator norm 

-.A 



a|||<C2||A|U^(,.(,^))^->^0. 



(8.27) 



(8.28) 



Moreover, being a holomorphic function of y we have that, for x large enough, (/ - C^^) is 
invertible and that its inverse is also a holomorphic function of 7. In particular, (8.25) imphes 



dyu^^(i-c^yo(d,ci^)[u^] 



(8.29) 
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A straightforward induction shows that there exist some coefficients 4 e Z such that 

''=1 ^=lPi=" 

° - C^)"' ° o---o(l- o [n,] . (8.30) 

This expression simpUfies at y = as A\^^ - and n+|^^g = I2. Hence, 

^"n.! =t Z ^''■'M^r^C^n)°(4^CtJo...o(5^^Cy[/2]. (8.31) 

We can slightly deform the different contours En, so as to regularize the explicit integral rep- 
resentation for the above chain of operators. Namely, recalling the construction of the jump 
contours for the matrices occurring in the different transformations appUed to the RHP for x, 
we are able to write 

^rn^l =Z Z c«^'''te'C^a))°fcc^a))°---°kc^(.)^^^^ (8-32) 

There, the contours l!^ are such that the - side of e|^j~^^ is at small but non vanishing distance 
from the + side of l!^, with the exception of a finite number of points of intersection (cf. 
Fig. 6). The matrix A corresponding to the contour l!^ is equal to M*^ - I2 on r^'\ to P - I2 

on dD-q Si^ and \o T - h on dDq Sr We emphasize that, already in (8.32), one can use the 
integral representation for the Cauchy operators without turning to boundary values. Indeed, 
the integrand appearing in (8.32) is already integrable on 1^'^^ = x ■ ■ ■ x 




p'O-i) 



Figure 6: Encased Contours xlt^ (in the case p = id). 



Finally, one infers from (8.32), from (8.22) and from the integral representation for 
(8.24) that, when A e r([-^;^]), there exist some recursively computable coefficients c^'''''^ 
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such that 



,1) r d'z 8('P"»(zr;x)...8(iP"»(zi;^) 



U - zif n fe-i - ze) 
e=2 



m=l 



J 



IPlm 



ze-i^ 



(8.33) 



In this expression, the integration is performed over the skeleton t^'^ = x ■ • • x xj:^ and 
the second summation is performed over integers ptj, I < £ < r, \ < j < n, with 1 < pm < n 
and < pcj < « for j > 1, and such that ^^f^iPt - n, in which we have introduced the notation 

pe - pio + ^=iiPij- 

• y-derivatives of 11"^ 

All the above observations also hold for the inverse matrix 11" ^ {X). Indeed satisfies the 
integral equation 



in which 



'Cl^[M} = C^^[VM], VMeM2(L'(2n)), 



(8.34) 



(8.35) 



and the matrix V is defined^ by the equation /2 + V - {I2 + A)~^ In other words, V is the 
adjugate of A (we remind that we consider 2x2 matrices and that det [/ + A] = 1). Hence, one 
easily sees that, for n > 1, 



p«^'-»(zi;x)...g«P"')(z?;;c) 
t 

i^-zi) n (zi-i - Ze) 

(=1 



X 



n 



F'" fe) f] 



m=l 



F>-{zt)-F'^i}i) 



pern 



(8.36) 



where q^'^P'*^^ {z; x) is the adjugate matrix of d^*^"'-* (z; x). 
• y-dcrivativcs of log det [/ + V]™*' 



^We stress that this matrix V has nothing to do with the differential operator usually denoted by V. 
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From the expressions (8.18), (8.33) and (8.36), it is easy to see that there exist some combi- 



natorial coefficients (£ 



({paY) 



; (with eJo"'^ = -<^pio.n) such that 



r>l, f>0 Tr+z r III -o .oil vv i 



r([-9;?]) 



/ dx' tr«f^'\{z;};x')[g] 



X 



r r+f 

iA-zir(A-Zr+i)Uiz{-i-zd n (z^-i-z^) 

^=2 ^=r+2 



{=1 



m=l 



1/ 

/ 



dfj. 



Pem\ 



(8.37) 



in which, in the terms t = 0, the empty product (A - Zr+\) YlfL'r+i ^^■f-'^ ~ ^''^ should be under- 
stood as 1. In this expression, tr[.'^'^'\{zj}; x')[g] corresponds to the following trace: 

tT%'%; x')[g] = tr jcTs ff;[Ao(z; x') - ^ A^^^^z; x')] | J , (8.38) 
ii?''\zj]; x')[g] = tr {6«^"»(z,; x')... 6«^"l\zi ; x') crs 

X g^^P^-"^\zr+i ; x') . . . ^>^'^'-«''\zr+?; x')} if r + ? > 1 . (8.39) 

In (8.38), Aq x') corresponds to the first term in the asymptotic expansion (5.1) of Ao(z). 

Remark 8.1. We have gathered in the term r - 1, ? = the contribution of the second term in 
(8.18), as well as the term that would correspond to the contribution of only one jump matrix A 
in the first term of (8.18). Note that, in this term r = l,t = 0, the only non-zero contribution 
come from the diagonal elements of A, hence from the sequence pio = n, pij = for j > 1 
(indeed we have gJq"'-* = -<5pio,n)- 

Remark 8.2. It is easy to see from these expressions that the integrals over x' are convergent. 
Indeed, it follows from the asymptotic expansion of the matrices Aq that 



5;;[Ao(z;x)-iAi^>(z;x)]|^^^=o(i^" 



(8.40) 



uniformly on the integration contour, so that an integration of the trace (8.38) is convergent. We 
emphasise that the trace (8.39) is at least O ((log xj^jx^) uniformly on the integration contour: 

indeed, each of the matrices S*^'^'^'^ or ^^'^o') is uniformly a 0((logx)^"' /x); in addition, the 
trace (8.39) involves a product of at least two such matrices since r + t > 2. These estimates 
guarantee that the integrals over x' in (8.37) are well defined. 
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8.2.3 Application of the density procedure and proof of Theorem 8.1 

In order to be able to apply the density procedure, we should express more explicitely the 
functional dependance of 5^ log det [/ + yC^'^)]^"'' I^^q on !ri*^'<*\ 

The F-dependence of 5" log det [/ + Vf^^ |y=o has already been explicitely extracted in 

(8.37), and all the g-dependence is contained in the traces tr|.'f'''\|zy}; x')[g]. Using the struc- 
ture of the matrices S^'''"'^ (z; x) and p^''''^''^ (z; x), one can be more precise concerning this g- 
dependance. Indeed, it follows from (8.23) that there exist some coefficients d|.|^^^'''*^^'\{z^}; x) 
which are piecewise smooth on the integration contour such that 



(8.41) 



Ee,=0 



Note that these coefficients D^.'^/^ ' ^' {{Zj}', x) are at least O ((log x)" /x^j uniformly on the inte- 
gration contour. Integrating these coefficients with respect to x, and defining 

X 
+00 

which are at least O ((log x)" /x) uniformly on the integration contours 1^^, one gets 



(8.42) 



5;jiogdet[/ + y]-''| = 2 Z ^ S^^^) \w^r 



r>l,f>0 t'' , r([ -?;?]) 



X 



E 



r r+t 
eu...,e.M±m {X - Zlf {X - Zr+l) D {Z{-1 - Zf) W {Zf-l - Z{) 
h £. =0 t=2 e=r+2 



X 



n 



m=l 



/ 



Z£-fi 



PCm 



(8.43) 



We stress that each e*^*^^*^ may only appear in combination with at least one F {z{) (as p(Q > 1): 
F (xt) ■ This guarantees that the functional above is continuous with respect to the sup 
norm on the space of symmetric functions in n variables z and n variables X. 

Before applying the density procedure, let us introduce one more useful notation. Define 
the finite difference operator 6^'"^(/i) by its action on pure product functions 

5^ V) • F\z) = F\z) - f''-'\z) FP(p) . 
This action naturally extends to symmetric functions of n variables 



(8.44) 



6f()U)-r„ 
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We remind here that pf - pfo + IP^^spfg, with pf - n. We have moreover used the 
notation [zeY' , which means that the variable zi is repeated pi times, and {fcl^'j. j , which 
means that the variable zi is repeated p\ times, zi is repeated p2 times, Zr is repeated 
Pr times. The purpose of introducing such tinite difference operator is to recast products of 
functions F"* {ze) - F'" (jj) appearing in (8.43) into a more compact form. Namely, 



n 

m=l 



-<l 



F'"izd-F'"(p) 
zt-p- 



dp 



P(m 



m=l 



/ 



dp 



Ptm 



n Pern k A n P(m 

m=l j=l i ^^''"'^ I m=l j=l 



■ FP'-P'' izd 



Therefore, setting e^'^^^ - <f>(z) and F(z) - (p(z) <f>(z), we get 

5;;iogdet[/ + y(^'^)]^"''l = 2] Z ^ T-P^^^ [ 



Pe 
Peo^ 

7C[{p/,).{^/)l 



D 



(1^ /1; -v) 



eu....er.M±m (A - ZlV U - Zr+l) U {Zt-l ' Zt) 11 - Ze) 

2i e;=0 {=r+2 



r+t n Ptm 

xnnn 

(=\ m=l 7=1 



f -^^^ ^'^'iHe,mA ■ n [P''(^^) r^'M . (8.46) 

J Zf- pe,m,j ^=1 



It follows immediately from the density procedure formulated in Theorem C.l that 7^"'' 
can be extended into a Unear functional on SymJJJ, W). Its action on a holomorphic function 
Tn e SymJJJ, W) is given as 



X 



r,x(r+<) 



ot)" 



e,,...,.,«el±l,0) (i - (i - Z.+ i) n (z^-1 - Z^) n {Z{-i - Zt) 
2 ei=0 ^=2 ^=r+2 



r+f Jl Ptm 

xnnn 

^=1 m=l ;=1 



/ 



Z{ - P{,m,j 



^^^^(Mf,m,j) 



{fer}i<^<.« 



. (8.47) 



The sum appearing in (8.47) is tinite, and since each integrand is a 0(log"x/x), I^^ [T^n] is 
itself a O (log" x/x). Hence Theorem 8.1 follows directly, since 



sub I 



(8.48) 
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8.2.4 Asymptotic expansion of I^^iT'n] and proof of Theorem 8.2 



In order to prove the existence of an asymptotic series of for I„ VFn\, i.e. for I^^{Tn\, we should 
be more precise on the structure of the coefficients (, i.e. show that they themselves admit an 
asymptotic expansion. Let us recall that these coefficients are obtained from the traces (8.41) 
involving the matrices {z(', x) and ^({''^j)) {z(', x). The latter being obtained from the jump 

matrix Aq. 

Clearly, all terms corresponding to an integration on the contours T^^ yield exponentially 
small corrections. Thus in what concerns the proof of an asymptotic expansion we can only 
focus on integrations along the contours dD±q^Sr We decompose the relevant contour [dDq^s U 
dD-q^sT^'^*'^ into sums of elementary skeletons dD^q^g^ = £^^a-i^,5i x • • • X dDa-^^,q,s^^,, where 
each (Ti takes values in {±1}: 



(8.49) 



o-q,0 



The matrices 6'^'''^'^ (z; x) (8.23) admit an asymptotic expansion into inverse powers of x on 
dDq s U dD-q s. This fact follows from the asymptotic expansion of Ao(z). The latter is obtained 
by taking adequate a-derivatives at a = or 1 of the asymptotics series (A.6) for ^(a, l;z) 
when z — > CX3. This is licit as, for fixed M, the 0(z~^"') estimate in the asymptotic series (A.6) 
is uniform with respect to a and since we perform a finite number of derivative with respect to 
a. This asymptotic expansion takes the following form: 



6«^^»(z;x) = . 



M Ad crj vxp+*g(z)yi [sjj^^^''*^-' {z; X+) 1 

^ ^ z ^ + 



M Ad [tep-+i;fe)]/2 [6*i'''^''*\z; X_) 1 

^ ] — ■ 1 ^ + 



logP" X 
logPo X 



z e dDq^s , 



z e dD-q^s , 



(8.50) 



the corrections being uniform on the contours. There the diagonal entries of the matrices 
b^l'''^''^\z:,X±) are some constants (i.e. x and z independent), whereas the off-diagonal ones 
are polynomials of degree pQ in the variable X± = log [±xipiz) - p^)]- 

An exactly similar structure holds for g^^PJ^^ (z; x) as it is adjunct to b^^Pj^^ (z; x). Hence, 
on the skeleton dDp-q^s - dDcriq,Si X • ■ • X dD^^^^^q^Sr+t^ the trace (8.41) can be expanded in the 
following form: 



M+l . N 



N=2 ki,...,kr+t=^ ei,...,er+te{±l,0) 
I.ki=N Sei=0 



{Pij), {ei} 
{(ri},{ki} 



({z,};logx) 



exp i Z fo) + ixpo-e] 1 n n \ 



r+t 

n ipiZi)-Po-e) 



(8.51) 
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There we have expUcitly factored out the dependence on the oscillating factor exp{ixLfe{Pa-f}. 
In this expression, the coefficients T)^ t are piecewise smooth on the integration contour and are 
polynomials of degree 2^ le^l peo in log x. 
We set, for > 2, 



l_iV^(()) 



{Ptj}, {^i) 



X 

({z,};log.) = <r''^ f 



{peMei) 
{(ri},{ki} 



({z,};log/) (8.52) 



when 2 eepa-f = 0, and 



k=N 



apt, 

r,t 



,1) 



[ptj], [ei] 
[a-iUki) 



(te};logx)e'^'2^'-*^' (8.53) 



otherwise. Note at this stage that, due to the constraint H = 0, there exist some integer mi^O 
such that e"^P<^«^ = e'^^ip*-?-). 

We then insert the result of integration into the expression for 7™'' rearrange the asymp- 
totic expansion into decreasing powers of x and separate the oscillating and non-oscillating 
parts. We obtain 



iV=l N=2-^ \X / 



(8.54) 



where 



N+l 



AN; nose) 



t(A'; osc) 



r,f o-i,...,crr+,=± ki,...,kr+t=l 

1^.1= I. I. t t 



{<Ti\,{ki} 



(iogx)[r„] 



r+< 



r,? ai,...,ar+i=± s=2 ki,...,kr+,=\ 



{Pfj}Aei} 

{(Tihiki] 



(8.55) 



(iogx)[r„]. 



(8.56) 



In these expressions. 



E ^ E E E 

r,t l<r+t<n Pm,-,Pen ei,...,er+(e{±l,0} 

{pei\,M r>\,t>i) piQ>\,U*J^p(=n Zei=0 



(8.57) 
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and the functional 3^, t is given by 



{pej}, {e;} 
{(TiUki) 



z 

r+t 



X 



{pej}, {et} 
{o-i],{ki] 



r([ 

{{zi} ; log x) 



dD 



^=2 
r+t n P(i 



(A - zif {A - zr+i) n ize-i - ze) xl {ze-i - zd '=i ^^^^^ " ^'''^ 

t=r+2 



n 



nun Ji;r^^.(/^^-i)K«i{{.^^^^^^^ 



l<^<r+r 
}l<^<r+t 



. (8.58) 



The coefficients i^^j being polynomials of degree YjI \^e\ Pto in log Jc, this ends the proof of 
Theorem 8.2 concerning the existence of the asymptotic expansion of cyclic integrals to any 
order in 1 jx. 

As it is presented, the form of this asymptotic expansion may look quite involved. Note 
however that the integrals over the contours 5£>^J'J'^ in (8.58) can be computed; they are ex- 
pressible in terms of partial derivatives of the function Tn at ±q (see Appendix D). It is proved 
in Appendix D that the non- oscillating term 4^;"°^*^) [:r„] of order can be expressed in terms 
of derivatives of ;F„ of total order not higher than A^, whereas the order of such derivatives does 
not exceed N -I'm the case of 4^' "^'^^ [!F„]. This property will be usefuU in [31], when we sum 
up the asymptotic behavior of a whole class of cycle integrals of the form (8.1) to obtain the 
asymptotic behavior of correlation functions. To perform this summation we use the knowledge 
of the number of partial derivatives applied to We finally point out that the integral over A 
produces derivatives of the function p{A) evaluated at ±q. 



9 More general kernels 

In the applications to quantum integrable models, one sometimes needs to use some modified 
versions of the GSK. 

Consider the operator I + Ve acting on [ ; ^ ] with kernel 

Vg (A,iJ.) = (^) F Ou) 6' (A) ff (ju) , (9.1) 

2m [6 {A) - eijj.)] 

where e± and F are defined as in (1.7). We assume in addition that ^ is a biholomorphism of U 
onto its image, that 6 {[-q ; q]) G Rand 9(1/ n 'H^) c 

Then the asymptotic behavior of logdet[/ -i- V^] when x — > cx) follows from Theorem 2.1. 
More precisely, we have the following corollary: 
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Corollary 9.1. Let Vg be as above. Then 

logdet[/ + ye] = 2 r d^y(/l)log'[e_U)] + J]log 

'J n J- 



-1 



where 



v{X) = log (1 + yF{A)) , %{A; q) = exp 
2m 



J 



e{A) - ein) 



+ o(l), (9.2) 



e{A) - ein) 



(9.3) 



and, as before, p± = [d^p{A)] U=±^, v± = v{±q). 

Proof — The change of variables 6iA) - ^ maps the kernel Vg on the one of the GSK 

> 2in (^ - 77) 

This kernel acts on [ 9 (-q) ;6(q)] which is, a priori, a non symmetric interval. However, it is 
enough to apply the transformation A A-{6{q) + 9 {-q)) /2 so as to recover the symmetry of 
the interval. Then, it remains to enforce the inverse transformations on the asymptotic formula 
for the Fredholm determinant of V. □ 

Let us write explicitly the asymptotics (9.2) in the case of the kernel 

2in sinh(/l - ju) 

as it plays a crucial role in the analysis of the asymptotic behavior of the two-point functions in 
the massless phase of the XXZ Heisenberg chain [31]. In this case, equation (9.2) reads 



<1 H 

logdet [/ + Vsh] = 2 d^v(.l)log'[e_(^)] + ^ J" did// 



v'{A)v(n)-v{A)V(n) 



log- 



[smh(2q) p'^x] 



d^ 




tanh {aq - A) 



+ o(l). (9.4) 



It is clear that the last equation can be used in order to obtain an analog of the asymptotic 
expansion for multiple integrals of the type (8.1) where the rational functions z- AaxQ replaced 
by the hyperbolic sinh(z - X). Namely, let 



r([ -?;?]) -« 



{^1 
iz] 



x(p(zj)-p(^j)) 



!=! sinh (zj - Aj) sinh {zj - Aj+i) 



(9.5) 
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Then under the conditions of Corollary 8. 1 one has the following asymptotic estimate 



Conclusion 

We have obtained in this article the leading asymptotic expansion of the Fredholm determinant 
of the GSK. As we have mentioned, our main motivation is to apply this result to the asymp- 
totic analysis of the correlation functions of quantum integrable models, using in particular the 
asymptotic study of multiple integrals performed in Section 8. This will be done in a forthcom- 
ing publication [31]. 

Another development is to extend the above analysis so as to handle truncated Wiener-Hopf 
operators with symbols having Fischer-Hartwig type discontinuities. The corresponding results 
are published in [33]. 

Let us also point out some unsolved problems. One of them concerns the derivation of the 
asymptotics of the Fredholm determinant of the GSK via the method based on its derivative 
over endpoint q. It would be important to obtain a complete justification of this method, since 
it is rather powerful and at the same time relatively simple. 

Another problem is to prove the conjecture on the Z - v periodicity for the asymptotic ex- 
pansion of the Fredhohn determinant. If this property does hold, then all oscillating corrections 
can be obtained from the non-oscillating ones via a simple shift of v by integer numbers. This 
could lead to a much simpler way to compute sub-leading corrections for such determinants. 
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A Some properties of confluent hypergeometric function 

For generic parameters (a, c) the Tricomi confluent hypergeometric function ^ {a, c; z) is one of 
the solutions to the differential equation 

zy" + (c-z)y'-ay^O. (A.l) 

It satisfies the properties: 

• Differentiation: 

r (a, c;z)^^[(a-c+ l)^(a + 1 , c; z) - Y(fl, c; z)] 

^^[(a-c + zma, c; z)-'P{a-\, c; z)] . (A.2) 

• Monodromy: 



Y(l-a,l;-z), (A.3) 



r2(a) 

where e - sgn(3(z)). In particular, 

^(a, Uze^"') = ^(fl, l;z)e-^"'" + — r-— Y(l - a, 1; -z), < 0, (A.4) 

^{a, 1 ; ze-^''^) = ^{a, 1 ; z)e^''"' - , ^(1 -a,\; -z), > 0. (A.5) 

• Asymptotic expansion: 

mc;z)~|;(-ir^^^^^^^z— z^oo, 4<arg(z)<^, (A.6) 

with {a)„ = T{a + n) IT (a). 

We have the following recombination between the Tricomi CHF T (a, c; z) and the Humbert 
CHFcD(a,c;z) 

O (a, c; z) = -f^e'^'-Y {a, c; z) + l^J^-(-<=)*^^ {c - a, c; -z) , (A.7) 
I {c - a) 1 (a) 



where e = sgn (3(z)), and 

0(a,c;z) = y ^i-. (A.8) 



n=0 
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Such a recombination formula allows to obtain the asymptotic expansion of the Humbert CHF: 

Y{c-a)\ z / ^ nl{-z) ^ ' 



r(a) n!z ^ ' 

One can estimate integrals involving a product of two CHF as below, either by using La- 
place-type integral representations for the functions 3> (a, c; z) and *F (a, c; z) or applying the 
method given in [33]. The latter uses Erdelyi's representation of Laplace transforms of products 
of CHF in terms of Lauricella function adjoint to some asymptotic expansion of Lauricella 
function. In any case, the result reads: 
+00 

J df{e-"V(a;/)- ij = -2/a, (A. 10) 





+ 00 



dt je-^'^T (a; t) + \ + ^ \ = 2ia - a[ilf(a) + i// (-a)] , (A.ll) 



and the Riemann integrability of the integrands is part of the conclusion. We remind the defini- 
tion of the functions T{a;t) and ip{a;t): 

if (y; = <t> (-V, 1 ; -it) (D (v, 1 ; it) , (A. 12) 

T(y;/) = -0(-y, l;-//)0(y, I; it) + {d^O){-v, l;-it)(D{v, I; it) 

+ O (-y, 1; -it) (d^^) (y, 1; //) . (A. 13) 



B Three preparatory Lemmas 



Here we prove three preparatory integration lemmas used in Section 6. 

Lemma B.l. Let 'R(u, t) be a function of two variables defined on I X R"*", where I is an open 
interval ofR. containing 0. Suppose that the partial applications u ^ Kiu, t) are {I) for all 
but finitely many t's and that t ^ Kiu, t) is Riemann integrable uniformly in u, i.e.: 

Vp > 0, VM > 0, Vmo e /, 3i/ > such that 



-1-00 

^ dt [d\'R{u,t) - ^^HiuQj)] 



M 



u e^-v + uq;v + UQ[r\ I,k e : 
Then for g e '^^ (f) 

6 +00 

^ xg{t)'R{t,xt) dt^giO) J 'R(0,t)dt + o(\) 

where the small o(l) is with respect to the successive limit x6 —> +00 and (5 — > 0. 



(B.l) 



(B.2) 
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Proof — One has 



5 xS 



J xig(t)m,xt)-g(0)'R(0,xt))dt^ JdyJ dy[g(y)'R(y,t)]dt 



(B.3) 



xy 



Consider a function 



g : (y, a,b)^ J dtdy [g (y) K {y, f)] 



(B.4) 



on the compact set [0;(5] x R+ x R+. g is clearly continuous on the interior and the uni- 
form Riemann-integrability of H (y, f) guarantees that it is continuous in an neighborhood of 
(*, +00, *), (*, *, +00) and (*, +00, +00). Hence, \g\ is bounded, say by B, as continuous function 
on a compact set. Thus, 



6 x5 



JdyJ dy[g(y)'R(y,t)]dt 



< 6B, 



xy 

which ends the proof of Lemma B . 1 . 

Lemma B.2. Let g e (J) for some open interval I containing 0, then 



(B.5) 



xdt 
xt 



-g(0) log ^(5 + 0(1), 



(B.6) 



where o (1) stands with respect to the successive limits x5 +00 and 6^0. 



Proof — We have 

6 S t 

g(0)log(x(5+l) + Jdf J dy^ 



r xdt 

J'^'^TTVt^ 




xt 







= g(0)log(5x + 



o(l) + j dyg' 



(y) log 



x6+l 
xy+l 



(B.7) 



(B.8) 



But, 



< sup \g'\ xi6- log (6 + l/x) /x)^0, 

[0;5] ^-"^ 



(B.9) 



which ends the proof of Lemma B. 2. 
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Lemma B.3. Let k be defined in terms ofv as in (2.2), and set 



Hiq) = - J dAdu 

-1 



A- 



+ ^eve log K(€q;q). 



Then, 



(B.IO) 



2v+— [logA-(^;^)J -2y_— [log K(-q;q)\ = — . 

aq q aq 



(B.ll) 



Proof — Using (2.2), one can express the derivative of H (q) as 



dHjq) 
dq 



= 2v(e^) rd/i^^^^^^^ + ;^6y(6^)^[log/.(6^;^)]. (B.12) 



-1 



Thus, proving (B.ll) amounts to estabUshing the equaUty 



V €v (eq) [log K (eq; q)] - — — ^ ^Y^y i^q) f dfi— 



ieq)-v'(M) 
eq-n 



(B.13) 



The latter follows from an integration by parts: 
Y_t ^ T" [log i^q'^ q)] 



dq 

= E 

e=± 



eVf < 



2q 



, V+-V- r v(n)-Ve-v'^{ji- eq) 
H z + e I dju— 



iH-eq) 



, V+-V- y_e - Ve + 2e^y' 

+ + 6 



2q 



-2eq 



v'Ou)-v; 



(v+ - v-f 



+ 



This ends the proof. 



C The density Theorem 



Theorem C.l. Let U, W be two open neighborhoods of\—q\q \ and let J^n "r , I holo- 
morphic fidnction on U" X W", symmetric separately in the n variables A and in the n variables 



66 



z. Then, for any compact subsets K ( resp. P)ofU( resp. W) there exists a sequence [(pp, <ppj 
in 'H(K)x9{ (P) such that 

^« ( {z} ) " Z fl ^"^'^ 'I'P ^^'^ uniformly on K" X P" . (C.l) 

' ' p=0 i=l 

Proof — Let K and P be as above. 

Let X - K" X P"/ ~, where the relation ~ is defined as follows: iA,z) ~ iA',z') if there 
exists a couple of permutations {o-,n) e (5„ x ®„ such that (A"',z'') - (A',z'), where A'^ stands 
for (Ao-(i), . . . , /lo-(n)). Since (5„ X (5„ is a discrete group, its action on K" x P" is by definition 
proper, i.e. "^L^K^X P" 

{(0-, ;r) e ®„ X (5„ : L"''" n L = 0} is discrete. (C.2) 

This ensures that X is a compact Hausdorff topological space. Moreover the space 'if {X, C) of 
continuous functions on X is canonically identified with the space of continuous functions on 
K" X P" that are symmetric in the first or the last n variables. 

Define the subspace Sof^ {X, C) as the subset of functions f^^^''^^ of the form 

5^M = n^w^fc), (C.3) 



;=1 



where {(p, 4>)e'H{K)x'H (P), and let S be the C*-algebra generated by 5 . We have that 5 and 
hence S separates points in X. Indeed, let (A, z) and (ji, y) be any two representatives in K" x P" 
of two distinct points in X. Thus 

• there exists Ai e K such that exactly p of the n coordinates of the «-tuple A are equal to 
Ai, whereas exactly q of the n coordinates of the n- tuple fi are equal to Ai, with p ^ q; 

• or there exists Zi e P such that exactly p of the n coordinates of the n-tuple z are equal to 
Zi, whereas exactly q of the n coordinates of the ?i-tuple y are equal to Zi, with p t q. 

The situation is similar in the case of the first n and last n variables, therefore we only treat the 
first case. By Lagrange interpolation there exist a polynomial Q such that, for any coordinate 
Ak of A and any coordinate p^ of yu satisfying Ak ^ Ai and pk ^ Ai, 

QiAk) = Q(pk) = l and Q{Ai) = 2. (C.4) 

The function 

' ' p=i 

separates the projections of (A,z) and on X. Thus »S is a C*-subalgebra of ^ (X; C) that 
separates points. It then follows by the Stone-Weierstrass theorem that S = '^^ {X; C). 

Let 'J^n be holomorphic on U" x W" and symmetric in the first and in the last n variables. 

o o o 

There exists compact sets <z U and P^ <zW such that K <z and P c P^. Here, stands 
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for the interior of K^. Thus the restriction of Tn^o K"x P" also belongs to (X^; C), with 
Xe^K'^x P^J ~, and therefore there exists (^p,'^p)p^-^ in {K^, C) x (P^; C) such that 

Tn ( 1 = Z fl (-^'^ ^^i^ uniformly on ^T," X P^. (C.6) 

In particular the sequence converges uniformly to l^n on (dK^)" X {dP^f, the latter set being 
compact. Therefore we have 

a/f, ap, n (Mi - A) (y< - Zi) 

uniformly in (A,z) e K"x P"^. Moreover, 

r and <pp(z)^ — (C.8) 

are holomorphic in K, resp. P. □ 



D Form of the sub-leading terms in 7^"*' 

In this appendix, we focus on the general structure of the sub-leading asymptotics of cyclic 
integrals. We show that the 1/x^ term in the non-oscillating part can be obtained as an action 
of at most partial derivatives of the function followed by an evaluation at ±q or by an 
integration over [-q; q]. 

In principle, the contour integrals defining (8.58) can be computed to the end. However, 
the result is quite intricate, and we do not need, for the further applications, the formula in its 
whole generality. Indeed, we are interested in a particular sub-class of such integrals. More 
precisely we shall focus on the sub-class that is susceptible to produce the highest possible 
derivatives of the function Here, by highest derivative we mean the total degree of all the 
partial derivatives that might act on the integrand. This subclass is identified in the upcoming 
lemma. 

Lemma D.l. Let r,t e N with r + t > I label negations cri,. .., (Xr+t e {±}. Also introduce 
sufficiently small numbers < 6i < • ■ ■ < Sr+t < q as well as positive integers ki,. . . ,kr+t. 
Finally, let G e 9i{pa-^q^Si X ■ • • X £)o-,+,9,<5,+,) and 
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Then, the integral [G] can be computed as some combinatorial sum involving 

derivatives ofG at the points aiq, the maximal order of such derivatives being equal to 

r+t 

^ ki -nr-n, + 6rfl + (5<,o - 2. (D.2) 
/=i 

Here nr, resp. nt, is the number of times the sequence (cri , . . . , (Tr), resp. (cTr+i, . . . , cr^+f), 
changes sign, and 6r,o, <5f,0 denote the usual Kronecker symbols. 

Proof — Let us prove the claim by induction on r + f. 
First, for r + ? = 1, (D.2) is obviously satisfied. Indeed, 

[G] = g'^f [G] = [d^-'G) (crq). (D.3) 

Let us now assume that the result holds for any function G up to some value of r + f. We 
will prove that it also holds for r + ? + 1. 

Note first that ^J7-'''*''^[G] - ^^Jf'' '^''^[G], in which G, {a-,}, {h} are obtained from G, {(t,), 
[ki} by a reordering of the variables. Hence, it is enough to prove the claim for ^r+i,t- We will 
have to distinguish two cases, depending on whether r+l = lorr+l>l. 

In the case r + 1 = 1, it is easy to see that 

which means that ^^^^f'^'''*^''^ [G] can be expressed in terms of derivatives of G of maximal order 

( H/=2 ki - \ -nt + 6tfi) + (h - 1), hence the result. 
Let us now consider the case r + 1 > 1 . We have 

MicTiUk}) ^ 1 / n — n n n 



o-g,o 



' 1^1-^2/ 



zi=criq 



f 



a-q,o 



n ^ n ^ n —L- n ^ 



At this point one should distinguish between the two possible cases: o-\(r2 = 1 or cricra = -1. 
We first assume o-\o-2 = 1 (i.e. that there is no change of sign between cri and cr2), and set 

Gl{Z2, ■ ■ ■,Zr+t+l) = ^ziG(zi, . . .,Zr+t+l) \zi=a-ig- Then 
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The latter can be expressed in terms of derivatives of G of maximal order k + {k2 + ki - k + 

ZLT' - -n, + 5,0 - 2) = Y^ili^' ki-1- nr^x -nt + 6,,o. 
We now assume that crjcra - -1. This leads to 



^1-1 



where the function 



*:=0 



k\ " 



Gu{z2,...,Zr+t+l) = 



d^^^Gizi,...,Zr+t+l) \zi-- 



CTiq 



(Z2 + Cr2^)*i *^ 



(D.7) 



(D.8) 



is holomorphic inside the integration contour dDa-2q,62 X • • • X dDtr,+,q,Sr+,- Once again, the result 
will be expressed in terms of derivatives of G and the maximal order of these derivatives will be 

^1 - 1 + ( Z;=r' -fir- fit + 6tfi - 2) = TjZI^^ ki ~ "'•+1 - + ^t,o - 2, which ends the proof 
of Lemma D.l. □ 

Remark D.l. The integral can be explicitely computed using the recurrence formulas (D.6) and 
(D.7). In particular, in the simplest case cri - ■ ■ ■ - (Xr and cTr+i - ■■■ - cTr+u we have 



r+t . 



(D.9) 



Ml,...,Mr+» {=1 



\Zi=(Tiq 



in which the parameters ue are summed over sets T( defined as 

r^ = lo,...,J]kj-J]^j-i\, r,,,= 0,..., 2 Z "^-1 

^ 7=1 M ^ ^ j=r+l j=r+l ^ 

r r r-\ X r r+t r+t-l \ 



(1 < ^ < r), 
(D.IO) 
(D.U) 



Corollary D.l. The subleadings terms of order N in the asymptotic expansion (8.54) /or the 
cycle integral InlTnl are obtained in terms of derivatives of the function Tn- More precisely, 

the non-oscillating term involves derivatives ofj^n of total order at most equal to 

N, whereas the oscillating one 1^^' involves derivatives ofTn of total order at most equal 

toN-2. 

Proof — In order to apply Lemma D. 1 to the integral over dD^^*^^ in (8.58), let us set 



(s) 



G({z))^ 



{{zi] ; log x) 



n 



ze - o-tq 



{X-Zxf{X-Zr+\) \}x\P^^^-P<^t 

nnn 



l=\ m=\ j=l 



/ 



W,m,j) 



Ze - pe,m,i 



ke 



[{Zcr],<t<r+t 



(D.12) 
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The poles at zi - A and Zr+i - A being outside of the skeleton dDj^^ ', this function is indeed 

holomorphic in a vicinity of the poly disc D^^^g\ 

Applying the result of Lemma D. 1 to this function and using the fact that Yjh = N + \ in 
(8.55), it follows immediately that the expression of I'i' '^°^'^\Tn\ cannot involve derivatives of 
the function Tn of order higher than A^. This maximal order of derivatives corresponds to ? = 
and V / (Ti = cr with cr = ± in (8.55). 

Similarly, as in (8.56) X^;' ^ ^» li^'"'^^\'^n] cannot involve derivatives of f„ of order 
higher than N - I. Moreover, due to the constraints H Q = and 2 QPa-f ^ 0, it follows that 
the variables cr,- have to take both values + and -, which means that either ? > 1 or ^ 1 in 
(D.2) (we recall that r > 1). Hence I^'°^''\'Fn] caimot involve derivatives of 'Fn of order higher 
than N -2. □ 
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